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ABSTRACT: This paper introduces a random statistical scan over the high-energy initial parameter 
space of the minimal SUSY B — L model-denoted as the B — L MSSM. Each initial set of points 
is renormalization group evolved to the electroweak scale-being subjected, sequentially, to the re¬ 
quirement of radiative B — L and electroweak symmetry breaking, the present experimental lower 
bounds on the B — L vector boson and sparticle masses, as well as the lightest neutral Higgs mass of 
~125 GeV. The subspace of initial parameters that satisfies all such constraints is presented, shown 
to be robust and to contain a wide range of different configurations of soft supersymmetry breaking 
masses. The low-energy predictions of each such “valid" point-such as the sparticle mass spectrum 
and, in particular, the LSP-are computed and then statistically analyzed over the full subspace of 
valid points. Finally, the amount of fine-tuning required is quantified and compared to the MSSM 
computed using an identical random scan. The B — L MSSM is shown to generically require less 
fine-tuninng. 
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I. INTRODUCTION 

The minimal supersymmetric standard model (MSSM) is the simplest possible N = 1 supersymmetric 
extension of the standard model of particle physics. It has the gauge group SU(3)c x SU( 2)l x 1/(1 jy, 
replaces each gauge field by a vector supermultiplet, each matter fermion by a chiral supermultiplet, has 
a conjugate pair of Higgs chiral superfields and has no right-handed neutrino multiplets. The MSSM was 
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introduced in various contexts in [1, 2] and has been extensively reviewed in [3]. However, without fur¬ 
ther modification, the most general superpotential for the MSSM contains cubic interactions that explicitly 
violate both lepton and baryon number and lead, among other things, to rapid proton decay-which is un¬ 
observed. The traditional solution to this problem is to demand that, in addition to the gauge group, the 
MSSM be invariant under an ad hoc Z 2 finite symmetry-/?-parity-which acts on individual component 
fields as (—i)3(B-L)+2s_ Although there have been many attempts to explain how /(-parity can arise in the 
MSSM [4-8] or be spontaneously broken [9-11]. These are all “non- mi nimai” in the sense, for example, 
that they require additional matter multiplets, have other ad hoc assumptions and so on. Is there a more 
natural and minimal approach to /(-parity in the MSSM? 

One begins by noting that, within the context of supersymmetry, /^-parity is a finite subgroup of 
1/(1) B _ L -see for example [12], It follows that /(-parity might arise as a consequence of demanding that 
the MSSM be invariant under a global U (1 )b~l symmetry. This would be consistent with current bounds 
on baryon and lepton number violation, and can also be imposed on the MSSM extended to include three 
families of right-handed neutrino chiral supermultiplets . However, one of the implications of the standard 
model is that a global continuous symmetry group is likely to appear - in its local form-that is, as a gauge 
symmetry. With this in mind, one might ask if /(-parity could arise as the consequence of extending the 
MSSM gauge group to include a new gauged 1/(1) b~l symmetry. It has long been known that the standard 
MSSM is anomalous with respect to gauged B — L symmetry, whereas the MSSM extended by three 
families of right-handed neutrino supermultiplets is anomaly free. Furthermore, this is the minimal such 
extension of the MSSM. We will call this anomaly free, minimal content theory the B — L MSSM, and 
propose that this is a more “natural” way in which /(-parity can arise in low energy supersymmetric particle 
physics. 

The B — L MSSM was identified from a low-energy “bottom-up” point of view in [13-15] 1 . In several 
papers [17, 18], these authors explored its structure and some phenomenological consequences. In addition, 
these ideas, and other directions, have recently been reviewed in [19]. Interestingly, the B — L MSSM was 
also discovered from a high-energy “top-down” viewpoint in [20-22], where it was shown that this model 
arises within the context of heterotic superstring theory [23, 24]. More specifically, the B — L MSSM is the 
low-energy effective theory associated with compactifying the Eg x Eg heterotic string on a Schoen Calabi- 
Yau threefold [25] with a specific class of SU( 4) vector bundles [26]. An important aspect of this high- 
energy point of view is that the parameters of the theory are specified near the gauge coupling unification 
scale, and then run down to the electroweak scale using the renormalization group (RG). This allows one 
to explore fundamental aspects of the theory-such as B — L and electroweak symmetry breaking. First 
1 See [ 16 ] for a similar idea in the context of Eg. 
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steps in this direction were taken in [27-29], where it was shown-for a restrictive set of initial parameters- 
that radiative breaking of both of these symmetries can indeed occur. A further study of the Wilson lines, 
spectra, mass scales and the unification of gauge couplings from the high-energy superstring point of view 
was presented in [30]. Combining the bottom-up and top-down approaches to the B — L MSSM, various 
aspects of both LHC and neutrino phenomenology were studied in the special case where a stop or a sbottom 
sparticle is the lightest supersymmetric particle (LSP) [31, 32]. 

However, all of the previous analyses involved specific assumptions-either about high-energy initial 
conditions or the low-energy structure of the theory. For the B — L MSSM to be a realistic contender for 
the low-energy theory of particle physics, it is essential that its initial parameter space be explored in a 
generic way, and that its low-energy predictions be compared with all present experimental data. This will 
be carried out in detail in this paper. We perform a statistical scan over a well-defined and wide range of 
initial parameters and, for each fixed set of such parameters, scale the theory down to low energy using 
the RGEs with specified threshold conditions. The results will be examined to determine the subset of the 
parameter space that, sequentially, 1) breaks B — L symmetry at a scale consistent with experiment, 2) 
breaks electroweak symmetry, 3) has all sparticle masses above their present experimental lower bounds 
and 4) predicts the mass of the lightest neutral Higgs scalar to be within 2a standard deviation from the 
ATLAS measured value of 125.36 GeV. A small subset of important results from this statistical scan were 
given in [33]. Here, we give present the details of the method use, as well as a wide array of new results 
and experimental predictions. 

This paper is organized as follows. In Section II, the B — L MSSM model at the TeV scale is described 
in detail. This includes a discussion of both B — L and electroweak symmetry breaking and much of 
our notation. The connection between the UV picture at the scale of gauge coupling unification and the 
TeV scale is outlined in Section III. Specifically, the important mass scales and energy regimes arc fully 
described and the relationships between the different mass scales arc presented. The content of these two 
sections is expanded upon in Section IV. Here, the details of the renormalization group equation (RGE) 
evolution of the parameters of the theory between the UV and TeV scales are discussed. This includes 
the appropriate input values of the parameters and the relevant equations-with reference to Appendix A 
when necessary. Special attention is given to the right-handed sneutrino RGE which drives radiative B — L 
symmetry breaking. The relationship between the different running parameters and scales-introduced in 
Section III—is further discussed. The latter part of this section outlines the experimental bounds used in our 
analysis. First, collider bounds arc discussed. These correspond to lower bounds on the physical sparticle 
masses-which arc closely related to the running mass parameters. The section finishes by describing-and 
implementing-the well-known bounds from flavor changing neutral currents and CP violation. 
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This paper approaches the connection between UV and TeV physics in a novel way. Specifically, instead 
of assuming some universal conditions or relationships between UV soft SUSY breaking parameters, we 
simply allow all such parameters to be within about an order of magnitude of some chosen SUSY breaking 
scale. Our approach will then be to scan all relevant SUSY breaking parameters, at the high scale, over 
this possible range. These parameters will then be RG evolved to the TeV scale, following the discussion 
in Section IV. This very general approach, as well as other details of our scan, arc described in Section V. 
Such an approach is especially applicable to the string realization of the B — L MSSM model, since, in 
that case, each chiral supermultiplet arises from a different 16 representation of 50(10). Therefore, their 
soft SUSY breaking masses do not obey boundary conditions at the high scale. Our analysis is also valid 
for a wide range of pure GUTs, which can impose disparate boundary conditions at the unification scale- 
including none at all. Finally a “meta” scan is conducted to choose the optimal value for the range of SUSY 
breaking mass scales. 

The optimal range, arrived at in Section V, is used to generate all subsequent results in this paper-the 
bulk of which arc presented in Section VI. Relating these results to the twenty-four phenomenologically 
relevant scanned parameters is daunting at best. Fortunately, a cohesive picture can be presented in terms of 
two so-called 5-parameters. These 5-parameters, which arc each the sum of the squares of SUSY breaking 
mass parameters, play an important role in radiative B — L symmetry breaking. Their role in this capacity is 
presented in Fig. 5. The results for all other experimental constraints can also be expressed in terms of these 
two parameters-see Fig. 6 and Fig. 7. A central result of this work is given is Fig. 8, which displays the 
frequency at which particles appeal - as LSPs in our scan. The section closes with histograms of the spectra, 
as well as some spectrum plots to help characterize specific features of our results. 

Fine-tuning in our approach is addressed in Section VII. While fine-tuning in the B — L MSSM is not 
drastically different than in the MSSM, there are several key issues to highlight. First, while one might 
expect that the scale associated with B — L symmetry breaking could introduce new contributions to fine- 
tuning, it is shown that this is not the case. Second, the MSSM, analyzed using the same methods as in this 
paper, typically yields equal or more fine-tuning than in the our model for similar initial points. Finally, an 
LSP analysis similar to Fig. 8, but with fine-tuning constrained to be better than one part in a thousand, is 
presented in Fig. 17. We conclude in Section VIII. 

In addition to the main sections of this paper, three Appendices are included to help elucidate various 
topics. Appendix A contains all one-loop RGEs for this model in the different regimes. The first part of 
Appendix B specifies how to relate the running soft mass parameters to the physical masses of the SUSY 
particles. The second part of Appendix B describes the procedure used to calculate the SM-like Higgs mass. 
Those readers interested in the details of how the random scan in this paper was conducted, are directed to 
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Appendix C. 


II. THE TEV SCALE MODEL 


Motivated by both phenomenological considerations and string theory, we analyze the minimal anomaly 
free extension of the MSSM with gauge group 


SU(3) C < 8 > SU(2) l 0 U{ 1) 3R 0 ■ 


( 1 ) 


As discussed in [30], we prefer to work with the Abelian factors U(\)zr ® U{!)b~l rather than U(1)y 0 
U(1)b- j- although they arc physically equivalent. This is motivated by the fact that the former is the 
unique choice that does not introduce kinetic mixing between the associated field strengths at any scale in 
their renormalization group equation (RGE) evolution. The gauge covariant derivative can be written as 

D = d — ihrigrtWri - i ———pbl-B' , (2) 

where / :j is the U( \ );»/,■ charge and the factor of 2 is introduced in the last term by a redefinition of the 
gauge coupling gnj- thus simplifying many equations. As discussed in [30] and throughout this paper, 
a radiatively induced vacuum expectation value (VEV) for a right-handed sneutrino will spontaneously 
break the Abelian factors U{1)sr x U(\)b~l to U{l)y, in analogy with the way that the Higgs fields 
break SU(2)l <S> U(1)y to U{1)em in the SM. For simplicity, we will refer to this as “B — L” symmetry 
breaking-even though it is technically the breaking of a linear combination of the U(1)^r and U(\)b~l 
generators, leaving the hypercharge group generated by 

Y = I.m + (3) 


invariant. The particle content of the minimal model is simply that of the MSSM plus three right-handed 
neutrino chiral multiplets. That is, three generations of matter superfields 


Q = 


L = 


u 

d 


v 

e 


1 (3,1,-1/2, —i) 

(3,2,0,-) ' ’ 

3 d c ~ (3,1,1/2,-2) 

( 1 , 1 ,- 1 / 2 , 1 ) 


( 1 , 2 , 0 ,- 1 ) 


along with two Higgs supermultiplets 

H u = 

H d = 


H+ 


H° 

u 


H d 

HT 


e c ~ ( 1 , 1 , 1/2,1) 

( 1 , 2 , 1 / 2 , 0 ) , 

( 1 , 2 ,- 1 / 2 , 0 ) . 


(4) 


(5) 
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We refer to this model throughout the remainder of this paper as the B — L MSSM. 

The superpotential of the B — L MSSM is given by 

W = Y u QH u u c - Y d QH d d c - Y e LH d e c + Y v LH u v c + nH u H d , (6) 

where flavor and gauge indices have been suppressed and the Yukawa couplings arc three-by-three matrices 
in flavor space. In principle, the Yukawa matrices arc arbitrary complex matrices. However, the observed 
smallness of the three CKM mixing angles and the CP-violating phase dictate that the quark Yukawa ma¬ 
trices be taken to be nearly diagonal and real. The lepton Yukawa coupling matrix can also be chosen to be 
diagonal and real. This is accomplished by moving the rotation angles and phases into the neutrino Yukawa 
couplings which, henceforth, must be complex matrices. Furthermore, the smallness of the first and second 
family fermion masses implies that all components of the up, down, and lepton Yukawa couplings-with the 
exception of the (3,3) components-can be neglected for the purposes of this paper. Similarly, the very light 
neutrino masses imply that the neutrino Yukawa couplings are sufficiently small so as to be neglected for 
the purposes of this paper. The //-parameter can be chosen to be real, but not necessarily positive, without 
loss of generality. The soft supersymmetry breaking Lagrangian is then given by 

-£ SO ft = Qm 3 s 2 + \m 2 w 2 + \m r wI + 1 -m bl b' 2 

+a u QH u u c - a d QH d d c - a e LH d e c + a u LH u D c + bH u H d + h.c)j 
+ iti 2 q\Q\ 2 + m\ c \u c \ 2 + m| c |d c | 2 + m||L| 2 + m| c |z> c | 2 + m%\e c \ 2 
+ m H u \H U \ 2 + m 2 H d \Hd \ 2 ■ 

The b parameter can be chosen to be real and positive without loss of generality. The gaugino soft masses 
can, in principle, be complex. This, however, could lead to CP-violating effects that arc not observed. 
Therefore, we proceed by assuming they all arc real. The a-parameters and scalar soft mass can, in general, 
be Hermitian matrices in family space. Again, however, this could lead to unobserved flavor and CP viola¬ 
tion. Therefore, we will assume they all arc diagonal and real. Furthermore, we assume that only the (3,3) 
components of the up, down, and lepton a-parameters arc significant and that the neutrino a parameters arc 
negligible. For more explanation of these assumptions, see Section IV B. 

Spontaneous breaking of B — L symmetry results from a right-handed sneutrino developing a non¬ 
vanishing VEV, since it carries the appropriate I^r and B — L charges. However, since sneutrinos arc 
singlets under the SU (3)c <S> SU (2 )l®U (l)y gauge group, it does not break any of the SM symmetries. 
To acquire a VEV, a right-handed sneutrino must develop a tachyonic mass. As discussed in [16, 34, 35], a 

Here and throughout this paper we use the term “tachyon” to describe a scalar particle whose m 2 parameter is negative. Al¬ 
though all m 2 parameters at high scale will be chosen positive, one or more can be driven negative at lower energy by radiative 
corrections. This signals dynamical instability at the origin-although a stable VEV may, or may not, develop. 



8 

VEV can only be generated in one linear combination of the right-handed sneutrinos. Furthermore, beyond 
the fact that its VEV breaks B — L symmetry, in which combination it occurs has no further observable 
effect. This is because there is no right-handed charged current to link the right-handed neutrinos to a 
corresponding right-handed charged lepton. Therefore, without loss of generality, one can assume that it 
is the third generation right-handed sneutrino that acquires a VEV. At a lower mass scale, electroweak 
symmetry is spontaneously broken by the neutral components of both the up and down Higgs multiplets 
acquiring non-zero VEV’s. In combination with the right-handed sneutrino VEV, this also induces a VEV 
in each of the three generations of left-handed sneutrinos. The notation for the relevant VEVs is 


(£ 3 ) = (4 = ~^ VLi ' ( H °) = ~^ Vu ’ ( H d) = 


(8) 


where i = 1, 2,3 is the generation index. 

The neutral gauge boson that becomes massive due to B — L symmetry breaking, Zr, has a mass at 
leading order, in the relevant limit that vr A > v, of 


M| » “ i (4 +s1l) 4 i 1 + sfkri) • 


(9) 


where 


V 2 = vj + v 2 u . (10) 

The second term in the parenthesis is a small effect due to mixing in the neutral gauge boson sector. The 
hypercharge gauge coupling is given by 

9 Y = 9 r sin Or = g B L cos Or, ( 11 ) 


where 


cos Or = 


9R 


( 12 ) 


\J9 2 r + 9rl 

Since the neutrino masses arc roughly proportional to the Y vi j and vj Jt parameters, it follows that 
Ypij < 1 and vj J% -C v u ^,vr. In this phenomenologically relevant limit, the minimization conditions 
of the potential arc simple and worthwhile to note. They arc 


4 = 


— 8 mb + g\ (yl - v 2 d ) 


9r + 9rl 


VLi =- 


y| (^i3 9 V d a *u i3 v u) 


m 2 f ~ iK-v 2 ) - 9 -f-v 


9bl„,2 


8 U R 


In o TO?r tan 2 B — m 2 rr 

-Ml =-r 2 + 


1 — tan 2 (3 


2b 


sin 2/3 


— 2 g + m H u + m H d 


(13) 

(14) 

(15) 

(16) 
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Here, the first two equations correspond to the sneutrino VEVs. The third and fourth equations arc of 
the same form as in the MSSM, but new B — L scale contribution to rri jj u and rriH d shift their values 
significantly compared to the MSSM. Eq. (13) can be used to re-express the Zr mass as 


m z r = -24 


1 + 


9r 


(17) 


9r + 9bL V Ry 

This makes it clear that, to leading order, the Zr mass is determined by the soft SUSY breaking mass of 
the third family right-handed sneutrino. The term proportional to v 2 /< is insignificant in comparison and, 
henceforth, neglected in our calculations. 

A direct consequence of generating a VEV for the third family sneutrino is the spontaneous breaking of 
/(-parity. The induced operators in the superpotential are 


7 ) €i L{ H u vn H d , 


(18) 


where 




7 f" aVR ■ 


(19) 


This general pattern of /(-parity violation is refeiTed to as bilinear i(-paiity breaking and has been discussed 
in many different contexts, especially in reference to neutrino masses- see references [36-39] for early 
works. In addition, the Lagrangian contains additional bilinear terms generated by vj Jt and vr from the 
super-covariant derivative. These are 


£ D " 2 VL * 


-Vr 


9 2 


2e i W + + v i W°) -g B LViB' 


-qrV^Wr + gBLV^B' 


+ h.c. 


( 20 ) 


The consequences of spontaneous /(-parity violation arc quite interesting, and have been discussed in a 
variety of papers. For LHC studies, see [17, 18] as well as recent work on stop and sbottom LSP’s in 
this context and the connection between their decays and the neutrino sector [31, 32]. Predictions for the 
neutrino sector were discussed in [16, 34, 35]. It was shown that the lightest left-handed, or active, neutrino 
is massless and that the model contains two right-handed neutrinos, referred to as sterile neutrinos, that arc 
lighter than the remaining two active neutrinos. Sterile neutrinos can influence the cosmological evolution 
of the universe due to their role as dark radiation. This effect was studied in [18]. 

In this section, we have focussed on the TeV scale manifestation of the B — L MSSM. However, the 
main content of this paper will be to study the connection between this low energy theory and its possible 
origins in Eg ® Eg heterotic suing theory, thereby linking some of the B — L MSSM phenomenology to 
high scale physics. In this context, our model is the remnant of an SO(\ 0) unified symmetry broken by two 
Wilson lines. Although the details of the various physical regimes of this theory, and the renormalization 
group scaling between them, were given in [30], we review them in the next section for completeness. 
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III. JOURNEY FROM THE UNIFICATION SCALE 

The goal of this section is to review the physics associated with the string construction of the B — L 
MSSM-from unification to the electroweak scale. After compactification to four-dimensions, the unified 
gauge group is 50(10). This is then further broken to the B — L MSSM gauge group by the turning on of 
two Abelian Wilson lines, denoted by X 3 R and Xb~l respectively. The energy scales associated with these 
Wilson lines need not be the same. In fact, exact gauge coupling unification at one-loop, which we will 
assume throughout this paper, requires that the scales be different- implying there is a two-step symmetry 
breaking process from 50(10) to the gauge group of the B — L MSSM. This leads to an intermediate 
regime between the two scales associated with the Wilson lines. The particle content and gauge group in 
this regime depends on which Wilson line turns on first. Defining the mass scales of X 3 R and Xb~l as 
M X3R and M Xb l respectively, we find the following two initial symmetry breaking patterns. 

• M Xb l > M X3R : 50(10) -> SU(3) C 0 SU(2) L 0 SU( 2) R 0 U(1) B -l, the “left-right” model 

• M X3R >M Xbl : 50(10) —>• SU(4:)c<S)SU(2)l^iU(1)sr, a modified version of the “Pati-Salam” 
model 


In each case, the subsequent turning on of the second Wilson line breaks the intermediate model to the 
B — L MSSM. 

Reference [30] studied these two cases and found that gauge coupling unification dictates that the Wilson 
line scales should be separated by less than an order of magnitude; that is, the intermediate regime is not 
very large. It follows that the TeV scale physics has little dependence on which of the above two models 
inhabit the intermediate regime. For simplicity, we will carry out our analysis under the assumption that 
it is the first of these symmetry breaking patterns that occurs. Hence, the intermediate regime contains the 
left-right model. We then make the identifications 

Mu = M Xb l , the scale of gauge coupling unification (21) 

Mi = M X3R , the intermediate scale (22) 


These scales will be further discussed below. 
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In the intermediate regime, the particle content of the left-right model consists of nine copies of the 
matter family 

/ 

d c 


Q~(3,2,1, J), Q c = 1 > ~ 

d \ u c 


(3,1, 2, —4 


(23) 


(1,2,1,-1), L c = | |~ (1,1, 2,1), 

u c 

two copies of a Higgs bi-doublet, which contains the MSSM Higgs fields, 

fti, U 2 ~ (1,2, 2,0) , 

and a pair of color triplets 


(24) 


(25) 


He ~ { 3,1,1,2), H c ~ (3,1,1, -2) . 


(26) 


Once the second Wilson line turns on, the extra particle content integrates out and one is left with exactly 
the spectrum of the B — L MSSM. 

At this point, it is important to make a quick note on notation for the B — L gauge coupling. Thus far, 
we have discussed the gauge parameter qbl, which couples to 4(2? — L) charge. As is well known, this 
gauge coupling has to be properly normalized so as to unify with the other gauge parameters. We use g' BL 
defined by 


9bl ~ 



(27) 


to denote the properly unifying coupling. The parameter g' BL couples to y | (B — L) charge and will appear 
in the RGEs. For quantities of physical interest, such as physical masses, gsL will be used. 

To fully understand the evolution of this model from unification to the electroweak scale, it should be 
noted that there arc five relevant mass scales of interest, two of which were mentioned briefly above. All 
five arc described in the following: 


• Mjj'- The unification mass and the scale of the first Wilson line. We assume that all gauge couplings 
unify at this scale. That is, g 3 = g 2 = g R = g' BL = gjj. 


• M\. The intermediate scale associated with the second Wilson line and the symmetry breaking 
SU(2) R —y [7(1 );>/-,> - that is, the right-handed isospin breaks into its third component. Since the 
gauge coupling of SU (2 ) R slightly above M\ is equal to the U (l): S /,> gauge coupling slightly below 
Mi, we use g R for both SU(2) R and U{1) 3R . All gauge couplings have trivial thresholds at this 
scale. 
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• Mb~l• The B — L scale is the mass at which the right-handed sneutrino VEV triggers U(1)zr 
U(1)b~l —> U(l)y. Physically, this corresponds to the mass of the neutral gauge boson Zr of the 
broken symmetry and, therefore, the scale of Zr decoupling. Specifically 


Mz r = Mr-L i 


(28) 


where Mz R depends on parameters evaluated at Mr-l~ see Eq. (9). Substituting Eq. (9) into this 
relation yields a transcendental equation that must be solved using iterative numerical methods to 
obtain the correct value for Mr~l- 


At this scale, we also evaluate the hypercharge gauge coupling using its relationship to the gauge 
parameters of B — L and the third component of right-handed isospin. This is given by 


9i 



g R sin 0 R 


/5 , 

y g 9bl cos 


Sr , 


(29) 


where 


COS Or = 


9R 


(30) 


\JfjR + %9rl 

Note that Eq. (29) is just a restatement of Eq. (11) with gauge couplings properly normalized for 
unification, including a rescaled hypercharge gauge coupling g\ defined by 

[5 


9i = 


9Y 


(31) 


• M$\jsy' The soft SUSY breaking scale. This is the scale at which all sparticles arc integrated out 
with the exception of the right-handed sneutrinos. The right-handed sneutrinos arc associated with 
B — L breaking and, therefore, are integrated out at the B — L scale. While there is obviously no 
single scale associated with the masses of all the SUSY partners, we use the scale of stop decoupling 
given by 


M susy = y/m^ m t - 2 . (32) 

This scale is useful because when the stops decouple, the parameter that controls electroweak sym¬ 
metry breaking, that is, the soft H u mass parameter, effectively stops running- see [40] for more 
details. Like the B — L scale, the SUSY scale must be determined using iterative numerical methods 
because the physical stop masses in Eq. (32) depend implicitly on the SUSY scale. 


• Mew^ The electroweak scale. This is the well-known scale associated with the Z and W gauge 
bosons of the SM. We will make the identification 


Mew = Mz- 


( 33 ) 
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For correct electroweak breaking, one must satisfy the conditions 


2b < 2 /r 2 + mjj u + m 2 Hd , 

(34) 

b 2 > if + M 2 Hd )(ii 2 + M 2 Hu ) . 

(35) 


The first constraint guarantees that the Higgs potential is bounded from below while the second 
indicates that the trivial vacuum is not stable. 


With the relevant mass scales appropriately defined, we can now discuss the physical regimes that exist 
in between them. To begin with, we will be interested in the evolution of the gauge couplings-since our 
assumption that they unify will help relate these disparate scales to each other. We present below, for each 
regime, the slope factors b a appealing in the gauge RGE’s 


d_ _i _ _ 
dt° a 2vr 


(36) 


where a indexes the associated gauge groups. Note that while Mjj > M\ 3 > M b ~l , T7susy> the hierarchy 
between the SUSY and B — L scales depends on the point chosen in the initial parameter space. Each of 
the two possibilities will be addressed below. 


• Mjj — M\. This regime is populated by the left-right model discussed above. In this interval, the b a 
factors arc 


& 3 = 10 , b 2 = 14, b R = 14, b B _ L = 19 . (37) 

We will refer to this scaling interval as the “left-right regime” and, when required, denote the associ¬ 
ated ^-coefficients by b]; ]i . 

• M\ — max ( MsusY > Mb-l)■ This regime is populated by the B — L MSSM model. The b a factors 
in this case are 


63 = —3, b-2 = 1, bn = 7. b B —L = 6 . (38) 

We will refer to this scaling interval as the “B-L MSSM regime” and, when required, denote its 
^-coefficients by 6 j^ L . 

The remaining two regimes depend on which of the following two cases occurs: M b ~l > MgusY-the 
“right-side-up” hierarchy-and M$ usy > Mb~l ~the “upside-down” hierarchy. 


right-side-up hierarchy: 
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• Mb~l — Msusy'- In this case B — L has been broken but SUSY is still a good symmetry, thereby 
giving an MSSM-like theory-that is, the MSSM plus two light right-handed neutrino chiral multi- 
plets. Another possible deviation might occur in the composition of the bino-more about this late. 
In general, however, this is the MSSM. Specifically, the gauge couplings in this regime evolve like 
the well-known MSSM gauge couplings with b a coefficients 

b 3 = -3, b 2 = 1, h . (39) 

5 

We refer to this interval as the “MSSM” regime and denote the associated 6 -coefficients by 6 ^ SSM . 


• Msusy — A Aw' In this regime, one simply has the SM with two sterile neutrinos. It has the well- 
known slope factors 

19 41 

6 3 = -7, 6 2 = - bl = (40) 

o 10 

We refer to this as the “SM” regime and denote the 6 -coefficients by 6 ® M . 
upside-down hierarchy: 

• Msusy — Mr-l- Now B — L remains a good symmetry below the average stop mass, where 
we effectively integrated out the SUSY partners. The resulting theory is simply a non-SUSY 
SU(3)c <8> SU(2)l <S> U(1) 3 r <8> U(1)b-l model, which also includes three generations of right- 
handed sneutrinos-the third of which acts as the B — L Higgs. The slope factors are 

63 =-7, 6 2 = y, 6 = 6 BL = f- (41) 

• Mb~l ~ Mew- Here, again, we have the SM with two sterile neutrinos and the slope factors given 
in Eq. (40). 


Given the above information, and the demand that all gauge couplings unify, we can solve for a given 
mass scale in terms of the others. First consider the unification mass-corresponding to the scale at which 
the four gauge couplings become equal to each other. Practically, it is derived as the energy-momenta at 
which g 3 = g ->. As is well-known, this will not be influenced by any scale that acts as a threshold for 
complete multiplets of a minimal group that unifies SU( 3) and SU (2)-for example, 5(7(5). The B — L 
and intermediate scales arc both such thresholds. The B — L scale is a threshold for singlets of 5(7(5), 


that is, the right-handed neutrinos, while M\ is a threshold for six new matter generations, a pair of Higgs 
doublets and their 5(7(5) color partners. All of these particles fit into the 1 , 5 , 5 and 10 of 5(7(5)- see 
Eqs. (23)-(26). Working through the algebra of setting g 3 (Mu) = .92 (Mu) yields 


M v 


27T(a 3 -a 2 ) 

e a 2 a 3 



b f M ) 



-fe| M +b| M -feM SSM ) 



(42) 
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where the superscripts on the slope factors indicate their regime of relevance and the cq take their experi¬ 
mental values at M z [41]: 


a 3 (M z ) = 0.118, a 2 (M z ) = 0.0337, an(M z ) = 0.0170 . 


(43) 


Inserting all the coefficients, the unification scale becomes 


Mjj ~ 2.186 x 10 


16 f Msxjsy \ 

V 1 GeV ) 


0.0417 


(GeV) . 


(44) 


Similarly, the intermediate scale can be solved for by setting gn(Mu) = 9bl(Mu) = < 73 (4%) and 
using the relationship between the gauge couplings of hypercharge , B — L and the third component of 
right-handed isospin given in Eq. (29). The intermediate scale is found to be 


Mi = 


« 2 ) 5 ( b SU_ b SM\ 5 ( b MSSM_ b SM +b SM_ b MSSM\ 

e “i “2 M z y 2 1 ; M si } s 2 y 2 i i j 


M, 


(36^+26“ -56^)' 


u 


5 ( 6 B L - i) L R ) +2 ( i ,LR L _ i ,BL i ) + 3 ( i ,LR_ i) : 


,BL I 

R I 


Substituting for Mjj using Eq. (44) gives 


Mi ~ 1.835 x 10 


17 f MsusY \ 

v 1 GeV J 


-0.486 


(GeV). 


(45) 


(46) 


These relationships arc displayed in Fig. 1. 


IV. THE FRAMEWORK 

The purpose of this paper is twofold: (i) to analyze how the high scale parameter space is associated 
with the TeV scale phenomenology of the B — L MSSM and (ii) to introduce a new way of analyzing this 
question. More specifically, (i) is an attempt to understand the kind of general UV conditions that yield 
valid, experimentally consistent points at the TeV scale-especially conditions that lead to B — L breaking 
and a viable B — L/electroweak hierarchy, a topic which has not been widely studied before. Such a study 
requires a framework for evolving parameters through the relevant regimes. The skeleton of this framework 
was discussed in the previous section. This section will continue this discussion and review the well-known 
experimental constraints on the SUSY parameters. Utilizing these constraints, the next section presents the 
motivation and the strategy for our scan. This addresses (ii). 

The approach to the RG evolution of the parameters is similar to other such work, with several deviations 
that will be highlighted below. The RGEs of interest arc calculated using reference [42] and arc presented in 
Appendix A. Gauge couplings and gaugino masses arc evolved up to the unification scale. The remaining 
parameters, Yukawa couplings, sfermion mass parameters and a-terms, arc only evaluated in the scaling 
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FIG. 1. The unification and intermediate scales as a function of the SUSY mass scale. These two scales have no 
dependence on the B — L scale. 

regimes below the intermediate scale. This is because in the string construction considered here, the scaling 
regime between the unification scale and M\ contains six additional copies of matter fields as well as an 
additional copy of Higgs fields.We note that each component field of a given generation of matter originates 
from a different 16 of SO (10). This is important and will be discussed later. Since these new Yukawa 
couplings are unknown, RG running them through this regime would not contribute to the predictability of 
this study. In practice, we implement these calculations piecewise starting with the analytically tractable 
equations first. These are the gauge couplings, gaugino mass parameters and the first and second generation 
sfermion mass parameters, as well as all sneutrino mass parameters. We then numerically calculate the 
evolution of the remaining parameters. 

As is traditional, we begin by inputting the experimentally determined parameters-that is, the gauge 
couplings and Yukawa couplings derived from fermion nrasses-at the electroweak scale. The initial values 
of the gauge couplings were given above in Eq. (43). For the purposes of this paper, the SM Yukawa 
couplings, which are three-by-three matrices in flavor space, can all be approximated to be zero except for 
the three-three elements which give mass to the third generation SM fermions. We use the initial conditions 

y t = 0.955, y b = 0.0174, y T = 0.0102. (47) 


For details on relating fermion masses to Yukawa couplings, see [45]. Here the lower case y represents 
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Yukawa couplings in the non-SUSY regime. These can be evolved to the SUSY scale, both in the right- 
side-up hierarchy, Eqs. (A5) - (A7), and up-side-down hierarchy, Eqs. (A 8 ) - (A10). At the SUSY scale, 
one has the non-trivial boundary conditions 


Vt(M, susy ) = ^(-^ susy ) sin/3 

2 / 6 , t(Msusy ) = Yb, T (MsusY) cos/3. (48) 


The boundary condition at the B — L scale is trivial. Above the B — L and SUSY scales, the Yukawa 
couplings arc only evolved up to the intermediate scale utilizing the RGEs in Eqs. (A14) - (A16). 

The gauge couplings in the various regimes were discussed in the previous section. With those solutions 
in hand, the RGE evolution of the gauginos can be easily derived. Gaugino masses are inputted at the 
unification scale and evolved down. Naively, one might expect gaugino mass unification. However, this is 
not always the case-as has been discussed in a number of contexts, see for example [43, 44]. Therefore, 
and to be as general as possible, we impose no relationship between the different gaugino masses at the 
unification scale. The general RGE for a gaugino mass parameter is 


Jt M - = 


baCX-aMa 
27T ’ 


(49) 


where a indexes the gauge groups. These equations can be solved analytically. For the gauginos associated 
with SU (3)c, SU ( 2 ) l, U(1)zr and (7(1) b~l one has 


i\,r _ M a{ M u) n m 

M a \t) — . ot a (t). 


(50) 


Oi a (Mu) 

The bino, however, is a treated somewhat differently for each of the two possible hierarchies between 
the B — L and SUSY scales. For the right-side-up hierarchy, at the Mb~l scale, we have three neutral 
fermions that mix: the third generation right-handed neutrino i/|, the B — L gaugino (blino) and the I 3 r 
gaugino (rino). This is a direct consequence of 77-parity violation in the B — L MSSM. As we will see, it 
is possible for a neutralino LSP mass eigenstate to have a significant /A component. The mixing between 
the third-family right-handed neutrino and the (7(1) gauginos is described in the (i/g, Wr, B') basis by the 
mass matrix 


^ 0 — cos 6rMz r sin 6rMz r ^ 

- cos 9rM Zr Mr 0 

\ sin 6rMz r 0 Mrl / 


(51) 


This mass matrix neglects mixing with the Higgsinos through the electroweak breaking Higgs VEV. This is a safe approximation 
since the lower bound on the Zr mass implies that the electroweak Higgs VEV will be negligible compared to the third-family 
right-handed sneutrino VEV. 



18 


Due to the RGEs monotonically pushing the values of Mr and Mrl down, they will typically be signif¬ 
icantly lighter than Mz R . It is, therefore, instructive to perturbatively diagonalize this mass matrix in the 
limit Mr, Mrl At zeroth order, the mass eigenstates are 


with masses 


B = Wr sin Or + & cos Or 


v / 2 (,/ 3 

* = ^ 


— Wr cos Or + B' sin Or) 
+ Wr COS Or - & sin Or) , 


(52) 

(53) 

(54) 


Mi = 0, m v o a = M Zr , m u c b = M Zr . (55) 

At first order, the effect of adding Mr and Mrl back into the mass matrix is to give the bino a mass of 

Mi = sin 2 OrMr + cos 2 0rM bl • (56) 


This shows that, in the right-side-up hierarchy, between the scales Mr~l and Msusy we have the gauge 
group and particle content of the MSSM plus two right-handed neutrino supermultiplets-that is, the two 
sneutrino generations that do not acquire a VEV. Below the B — L scale, the bino mass is 


Mi{t) 


M\{M B - L ) 

(Mr-l) 


ai(t). 


(57) 


In the upside-down case, all neutralinos arc diagonalized at the SUSY mass scale. 

The running of the tri-linear o-terms is straightforward. Their initial values arc randomly generated at 
the intermediate scale. The a-term RGEs in the B — L MSSM regime arc given in Eqs. (A17) - (A19), while 
those for the MSSM arc in Eqs. (A20) - (A22). All relevant threshold conditions arc trivial. 

The RGEs for the square of the soft sfermion mass parameters can be broken into two categories: 1) 
those with simple analytic solutions-given in Eqs. (A29) - (A34) and Eqs. (A35) - (A40) for the B — L 
MSSM and MSSM regimes respectively-and 2) those requiring numerical solutions-given in Eqs. (A42) - 
(A48) and Eqs. (A49) - (A55) for the B — L MSSM and MSSM regimes. These parameters arc all inputted 
at the intermediate scale. The third generation right-handed sneutrino is then evolved to the B — L scale- 
while all other sfermion mass squared parameters arc RG evolved to the SUSY scale. The third generation 
right-handed sneutrino mass squared plays an important role here since, when it runs negative, it triggers 
B — L breaking as was discussed in detail in [27, 28]. The right-handed sneutrino mass RGE is 

167r 2 —m?c = —3 g 2 R L Mg L - 2gj i M^ i + -QrlSrl - 9rSr , (58) 


At some points in parameter space, it is possible that the required limit will not be satisfied and there will not be a mass eigenstate 
that can clearly be identified as the bino. However, since the scaling regime between Mb-l and M$usy a l wa Y s small, the 
errors introduced by assuming the existence of a bino are insignificant. 
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where 


Sbl = Tr (2 m 2 Q — m~ c — m| c - 2m| + m| c + m~ c ) , 

„ o 2 ^ / 3 2 3 2 1 2 1 

Sr = mfr — rnir + Tr —m~ c + - rri~ - m}, c + -m 

n u n r i \ o^ o 1 


+ 


+ 


1 9r( m i) 9r(Mb-l) , , ,,, A 2 , 9bl( m b-lj „ /f __, 


Mb(A %) 2 + r LV ; Mbl(MuY (61) 


(59) 

2 ~u~ • 2 "' , d c 2 ' 2' " e ^ 1 

Despite the lack of a large Yukawa coupling, the right-handed sneutrino mass can still be driven tachy- 
onic by appropriate signs and magnitudes of the S'-terms defined in Eqns (59, 60). To emphasize this, the 
analytic solution to the sneutrino mass RGE is presented here. It is 

ml 3 (M B - L ) = ml 3 (Mi) 

( A (H/TA ^ „4 
14 gfr 

1 9r( m i) - 9r( M B-L 

14 4(M0 16 gl^A'h) 

Recall that the value of any Abelian gauge couplings grows larger at higher scale. Therefore, we see 
that a tachyonic sneutrino is only possible when Sr{M\) is negative and/or S B l(M i) is positive. This 
demonstrates the central role played by the S'-terms in the breaking of B — L symmetry. Note that in 
typical unification scenarios all soft masses are “universal” and, hence, both S-terms vanish. However, it 
was mentioned earlier that, in this string construction, different elements of a given generation arise from 
different 16 representations of SO(IO). Therefore, the soft masses of a given generation arc generically 
non-degenerate. Hence, the S-terms can be non-zero. 

As mentioned above, Mz R — \/2|mm| and the relationship 


-Sfl(Mi) - 


1 9%l(M i ) g 2 BL (MB~L) 


Sbl(M{) 


Mz r {Mb~l) = Mr-l 


(62) 


is used to iteratively solve for the B — L scale. The SUSY mass scale must also be solved for iteratively 
using the equation 


\f m tT ( m s\jsy) irif 2 (M S usy) = M SV sy , 


(63) 


where mg < mg arc the physical stop masses. The relationships between the soft mass parameters and 
the physical masses arc given in Appendix B 1. The soft mass squared parameter for the up-type Higgs is 
driven tachyonic, as usual, by the large top Yukawa coupling. Furthermore, the decoupled values of the soft 
Higgs mass squared parameters arc used to calculate the //- and 6-terms using Eqs. (15) and (16). 

The soft mass parameters have non-trivial boundary conditions at the B — L scale due to the effects of 
the B — L and I%r D-terms: 

m l( M B- L ) ~ = -^(g 2 R + g 2 BL ) v} { (l 3R - Y sin 2 Or) 

- ~ m z r (hR - Y sin 2 d R ) , 


(64) 
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where Mg_ L and Mg_ L indicate a scale slightly below and slightly above the B — L scale respectively, 
and /;j /,< and Y arc the third component of right-handed isospin and hypercharge of a generic scalar <i>. 

Once the mass parameters have been properly evolved to their appropriate scales, the physical masses 
can be evaluated. For much of the spectrum, this has been discussed in the literature, see for example [3], 
and has been included in Appendix B 1. The new element here is the mass of the scalar associated with the 
third generation right-handed sneutrino-the B — L Fliggs-degenerate with the Zr mass. In addition, the 
calculation for the SM-like Higgs mass is crucial since the experimentally measured value of ~125 GeV 
requires substantial radiative corrections from the stop sector in the MSSM. In this paper we follow the 
approach of references [46—48 ]—tak i ng into account the decoupling scales of the two stops, matching the 
quartic Higgs coupling at those scales and RGE evolving the quartic coupling to the electroweak scale to 
calculate the Higgs mass. Full details are given in Appendix B 2. 

Once a given physical mass is calculated, it is compared to current lower bounds or, in the case of the 
SM Higgs, the experimentally measured value. If a given initial set of parameters predicts a physical mass 
that is inconsistent with current bounds, it is rejected as being an invalid point. These lower bounds are 
discussed in the next subsection. 


A. Collider Constraints 

The bounds placed by collider data on SUSY masses are, in general, model dependent. That is, they 
depend on the spectrum and decay modes. Despite the much larger energy of the LHC, LEP 2 still has 
competitive bounds on colorless particles that couple to the Z and/or the photon-including sleptons in 
scenarios with both //-parity conservation [49, 50] and violation [51], bounds on charginos [52, 53] and 
bounds on sneutrinos in the case of //-parity violation [51]. As one may expect, due to the relatively clean 
environment at LEP, these bounds are close to one half the center of mass energy of LEP 2. Therefore, for 
simplicity, we proceed with the bound that all colorless fields that couple to the photon must be heavier than 
100 GeV. That is, 

mj, m-± > 100 GeV, (65) 

where i is any charged slepton. Colorless states that couple to the Z, the left-handed sneutrino, must be 
heavier than half the Z mass: 

rriu L > 45.6 GeV, (66) 

Colorless states that do not couple to the Z, such as right-handed sneutrinos/neutrinos and the bino, have 
such small collider production cross-sections that they do not have collider-based lower bounds. Wino and 
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Higgsino neutralinos arc degenerate with their chargino partner, thereby effectively putting a lower bound 
of 100 GeV on those states as well. 

The bounds from the LHC are much more dependent on the parameters. For example, if one investigated 
the bound on, for example, degenerate squarks in this model with a neutralino LSP, those bounds could be 
significantly different than in the case of a sneutrino, or some other, LSP. Allowing the squark masses to 
split would further alter the lower bounds. In fact, a full treatment would involve calculating the signatures 
of a given point in parameter space, comparing the number of events to the most recent LHC bounds on 
such events, and determining if the parameter point is valid. We do not expect the details of these lower 
bounds to heavily affect our results. We will, therefore, simply use the naive bounds 

mg > 1000 GeV, m- g > 1300 GeV , (67) 

which are based on recent CMS [54] and ATLAS [55] studies of the /(-parity conserving MSSM. In these 
studies, the colored states decay into jets and missing energy-possible final states in our model whenever 
the LSP decays into neutrinos. In this paper, we impose these bounds except in the case of a stop or sbottom 
LSP. These two cases were explicitly studied in [31, 32] and yielded the following lower bounds: 

admixture (right-handed) stop LSP: rn^ > 450 (400) GeV, > 500 GeV, (68) 

where t\ (b\) denotes the lightest stop (sbottom). Here, right-handed refers to a stop that is almost com¬ 
pletely right-handed-that is, a stop mixing angle, 0 f - > 85° or, equivalently, a state composed of 99% right- 
handed stop-while admixture stop refers to all other stops. This distinction is based on the phenomenology 
of the stops; right-handed stops have significant decays into a top quark and neutrinos while admixture stops 
decay almost exclusively to a bottom quark and a charged lepton. 

The lower bound on the Zr mass from LHC searches is 2.5 TeV [56, 57]. Finally, we require that the 
Higgs mass be within the 2n allowed range from the value measured at the ATLAS experiment at the LHC. 
We naively obtain the two sigma range by adding in quadrature the systematic and statistical uncertainties 
from [58], and multiplying the result by two: 

m h o = 125.36 ± 0.82 GeV. (69) 

See [59] for comparable data from CMS. A summary of the collider bounds mentioned above is given in 
Table I. 


B. Constraints from Flavor and CP-Violation 

A large number of low-energy experiments exist which place constraints on the SUSY parameter space. 
Some of the oldest and most well-known are the constraints placed on flavor changing neutral currents 
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Particle(s) 

Lower Bound 

Left-handed sneutrinos 

45.6 GeV 

Charginos, sleptons 

100 GeV 

Squarks, except for stop or sbottom LSP’s 

1000 GeV 

Stop LSP (admixture) 

450 GeV 

Stop LSP (right-handed) 

400 GeV 

Sbottom LSP 

500 GeV 

Gluino 

1300GeV 

Zr 

2500 GeV 


TABLE I. The different types of SUSY particles and the lower bounds implemented in this paper. 


from the analyses in references [60-62]-for example, those arising from K — K oscillation-and on CP 
violation [63-67]-for example, from electric dipole moment measurements. Genetically, the implication of 
these constraints are, approximately, as follows: 

• Soft sfermion mass matrices are diagonal. 

• The first two generations of squarks are degenerate in mass. 

• The trilinear o-terms are diagonal. 


• The gaugino masses and trilinear a-terms are real. 


In addition, it is typically assumed that the soft trilinear a-terms are proportional to the Yukawa couplings, 
that is, generically a = YA for each fermions species. Each A is a dimensionful real number on the order 
of a TeV, while each Y factor is a dimensionless matrix in flavor space. This condition effectively makes 
all non-third generation trilinear terms insignificant. Note that this assumption does not immediately follow 
from the above experimental constraints. However, significant radiative corrections to fermion masses, 
proportional to the a-term, can arise in SUSY, as first discussed for fermions in references [68-70]. For 
example, a down quark mass is modified by gluino exchange, through the diagram in Fig. 2, as follows: 

AM d = M d Mg ^ ^ + utanfi) I , (70) 


where 


xy In | 
(x — 


+ yz In | + xz In 

y)(y - z){x - z) 


I(x,y,z) 


X 


(71) 
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and Mg is the gluino mass. If u,i is on the order of a TeV, this radiative correction can be quite large, possibly 
larger than the down quark mass. If this were the case, the radiative correction would have to be fine-tuned 
against the tree-level contribution to reproduce the correct down quark mass. This motivates allowing only 
the third generation a-terms to be significant. Hence, our assumption that 

o(Afi) = Y(Mi)A(Mi) . (72) 

This makes all a-paramctcrs, except for those associated with t, b and r, insignificant. For simplicity, we 
will choose all other a-parameters to be zero. 


m d ( + /ztan/3 

v * d w 


dr, / 


' d 


R 


9 


dL 


dR 


FIG. 2. Radiative contribution of the gluino to the down quark mass. Similar contributions exist for the other fermions. 

Summarizing the above, we employ the following constraints motivated by low-energy physics: 

= diag (m|, m? , m|) , q = Q,U, D , 
mj = diag , t = L, E , v c , (73) 

a f (M 1 ) = Y f (M 1 )A f (M 1 ) , / = t, b, r . 

Note that all of these constraints can be implemented at the high scale, since RG evolution to the SUSY 
scale will not spoil these relations. Furthermore, we do not assume here that the first and second generation 
slepton masses arc degenerate-unlike the squark masses- since this is not required by low energy experi¬ 
ments. The degeneracy or non-degeneracy of these states will not, however, greatly effect the results of this 
paper. 


V. THE PARAMETER SPACE AND SCAN 

The previous section reviewed the framework used in this paper for connecting the high scale to LHC 
accessible physics. It remains at this point to discuss the input values for the SUSY breaking parameters. 
In this section, we introducing a novel way to analyze the initial parameter space of a SUSY model. While 
there have been many studies of specific, fixed boundary conditions at the high scale, and some recent 
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interesting discussions of random parameter scans at the TeV scale [71, 72], a study that combines both 
has not-until now- been undertaken. Specifically, our approach in this paper is to make a statistical scan 
of input parameters at the high scale-followed by a RG evolution of each set of those parameters to the 
TeV scale and an analysis of which of these high scale initial conditions lead to realistic physics. While 
the soft SUSY breaking sector contains over 100 dimensionful parameters, the constraints of low energy 
experiments discussed in the previous section only allow collider significant values for about a fifth of these- 
24 to be specific. These, along with tan f3 and a discussion of the sign of certain parameters, are presented 
in Table II. 

The high scale initial values of the 24 relevant SUSY breaking parameters arc determined as follows. 
To conduct our scan, we make the assumption that there is only one overall scale associated with SUSY 
breaking. This assumption does not require that the soft mass parameters be equal to each other, or even 
have similar values. It does, however, require that these parameters be at least within an order of magnitude, 
or so, of each other. To quantify this, we demand that any dimension one soft SUSY breaking parameter be 
chosen at random within the range 

M 

( j,Mf ), (74) 


where M is the mass setting the scale of SUSY breaking and / is a dimensionless number satisfying 
1 < / < 10. We will further insist that any such parameter be evenly scattered around M; that is, that M 
be the average of the randomly generated values. Clearly, this will not be the case if parameters arc chosen 
from a uniform probability distribution in the range (y-, Mf )-referred to as a “flat prior”. Instead, a “log 
prior” is adopted. This means that the natural logarithm of a given soft SUSY breaking parameter is chosen 
from a uniform distribution in the range 


(In 


M 1 


1 GeV / 


■7 > ln 


M 


1 GeV" 


/ )• 


(75) 


With a log prior distribution, M is the geometric mean of the randomly generated parameters. In addition 
to the dimensionful soft masses, we also scan tan /3 as a flat prior in the range (1.2, 65), thus selected so that 
all Yukawa couplings remain perturbative through the entire range. Furthermore, we randomly generate the 
signs of //, the three tri-scalar couplings au,. T , and the four gaugino masses 

In this paper, we arc interested in the low energy spectra being accessible at the LHC or a next generation 
collider. Therefore, in addition to the experimental constraints mentioned in the previous section, we further 
demand that all sparticle masses be lighter than 10 TeV. We call any point that satisfies this, as well as all 
previous criteria, a “valid accessible” point. The parameters M and / arc chosen in such a way so as to 
maximize the number of such points. To determine the values of M and / which yield the greatest number 
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of valid accessible points, we begin by making a ten by ten grid in the M — f plane. At each of these 
hundred points, we randomly generate one hundred thousand initial points in the 24-dimensional parameter 
space discussed above, RG scale them to low energy, and count the subset that satisfies the experimental 
checks discussed above. We then plot curves corresponding to a constant number of valid accessible points 
in Fig. 3. The plot shows a broad peak or plateau, the center of which maximizes the number of such points. 
This maximum occurs approximately for 

M = 2700 GeV, / = 3.3 . (76) 

These values will be used to generate the results in the remainder of this paper. Note that for these values, 
the smallest soft parameter is maximally about an order of magnitude away from the largest soft parameter. 
Specifically, for M and / in Eq. (76), the ranges for the random scan of each parameter are given in Table II. 



FIG. 3. A contour plot in the M - f plane of the number of valid accessible points; that is, points that meet all 
experimental constraints given in the previous section and for which all sparticles are below 10 TeV. A broad peak or 
plateau is evident around M = 2700 GeV and / = 3.3. 

The existence of a peak in Fig. 3 around moderate values of / is a consequence of combining the 
various experimental checks we apply to each of the randomly generated points. For a fixed value of M, 
some individual checks favor higher values of /, while others favor lower values. This is analyzed in 
terms of the “survival rate”. The survival rate for a given check is defined as the number of points in the 
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Parameter 

Range 

Prior 

niq^ = mq 2 , m.q 3 : q = Q , u c , d c 

(820, 8900) GeV 

log 

m h > m t 2 > m l 3 '■ 1= L,e c ,u c 

(820, 8900) GeV 

log 

m Hu , m Hd 

(820, 8900) GeV 

log 

AH : f = t,b,r 

(820, 8900) GeV 

log 

\M a \ : a = R, BL, 2,3 

(820, 8900) GeV 

log 

tan /3 

(1.2, 65) 

flat 

Sign of n,a,f,M a : f = t,b,T a = R,BL, 2,3 

[-.+] 

flat 


TABLE II. The parameters and their ranges scanned in this study, as well as the type of prior. The ranges for the soft 
SUSY breaking parameters are optimized to produce the greatest number of valid points with all masses below 10 
TeV. 

24-dimensional initial parameter space surviving that check as a percentage of the number of points that 
survived all previous checks. This will be discussed in detail in the next section for the fixed values of M 
and / given in Eq. (76). Here, for M = 2700 GeV, we analyze the impact of the parameter / on the various 
survival rates. The peak around moderate values of / shown Fig. 3 can be understood by observing how 
the survival rates for different checks depend on /. This is shown in Fig. 4. The B — L symmetry breaking 
check and the Zr lower bound check both prefer higher values of /. This is because such values favor 
larger 5-terms and thereby promote B — L symmetry breaking. The EW symmetry breaking check favors 
lower values of /. Intuitively, this is not surprising since universal boundary conditions (which correspond 
to / = 1) in the MSSM allow electroweak symmetry breaking. The sparticle lower bounds check favors 
low /. This is because larger / leads to larger 5-terms which, in turn, can drive some sparticles masses to 
be light through the RGE’s. The Higgs mass check also favors low / because larger 5-terms may drive the 
stop masses away from the ~ TeV value favored by the Higgs mass. With some checks favoring large / 
and others small /, it is not surprising that all checks taken together favor a moderate values of /. 


VI. RESULTS 

All of the following results arise from a scan consisting of ten million randomly generated points with 
M = 2700 GeV and / = 3.3. We will refer to this as the “main scan”. Recall that in the previous section, 
in addition to the experimental constraints, we imposed an extra condition that all masses be lighter than 
10 TeV. Technically, this was done to ensure the maximum number of valid accessible points-that is, valid 
points with masses accessible to the LHC or a next generation collider. Having done this, we, henceforth. 
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FIG. 4. Survival rates of the various checks as a function of / for M = 2700 GeV. The B — L breaking and Zr lower 
bound checks favor larger / while the others favor small /. All of the checks taken together favor a moderate value 
of / ~ 3.3. 

remove this additional condition. That is, the only constraints in the main scan are the experimental ones 
given in Section IV. All valid points must satisfy these constraints as well as some other checks, which are 
reviewed here briefly. 

In order to be valid, a point must break B — L symmetry, the Zr mass must be above the lower bound, 
electroweak symmetry must be broken and since the stops play a crucial role in electroweak breaking, we 
designate a check that the stops are not tachyonic. Since our numerical analysis uses an iterative process to 
solve for the SUSY and B — L scales, it is possible that a point may pass a check on the first iteration but 
fail it on the final iteration, although this is very uncommon. Therefore, we include a “spill” check of the 
Z[{ bound, electroweak breaking check and non-tachyonic stop checks. Points that pass these spill checks 
did so on the final iteration of solving for the B — L and SUSY scales. Furthermore, a valid point must have 
B — L and SUSY scales that converge to a value in the iterative solution process. We also check that-in 
addition to the stops-all other SUSY sparticles are not tachyonic and satisfy the imposed mass bounds. 
Finally, we check that the Higgs mass matches its experimental value. 

All of these conditions are listed in the first column of Table III. The second column lists the number of 
points in the main scan that passed that check, out of ten million. The third column is the same information 
listed as a percent of the number of points in the main scan. The fourth column is the same information 
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listed as a percent of the number of points that passed the previous checks. We refer to this quantity as the 
rate of survival for each check. This is an interesting quantity because it quantities how easy or hard it is for 
a randomly generated point to pass that specific check. 


check 

number surviving 

percent surviving 

rate of survival 

B — L breaking 

2,225,704 

22.3 % 

22.3 % 

Zr bound 

919,117 

9.19 % 

41.3 % 

EW breaking 

722,750 

7.23 % 

78.6 % 

non-tachyonic stops 

619,668 

6.2% 

85.7 % 

Zr bound spill 

597,988 

5.98 % 

96.5 % 

EW breaking spill 

565,272 

5.65 % 

94.5 % 

non-tachyonic stops spill 

553,592 

5.54 % 

97.9 % 

convergence 

553,150 

5.53 % 

99.9 % 

sparticle bounds 

276,676 

2.77 % 

50% 

Higgs mass 

58,096 

0.581 % 

21 % 


TABLE III. This table shows all of the checks applied to the randomly generated points. It specifies the number of 
such points passing each check, as well as their percent of survival. The fourth column is the most informative because 
it provides insight into how likely it is that an individual check is satisfied by a randomly generated point. Because 
the SUSY and B — L scales are solved for iteratively, it is possible to pass a check in the first iteration and fail it later. 
A passed “spill” check indicates that that check was passed in the final iteration. 

A striking feature of Table III is that B — L breaking happens robustly. This was one of the central 
questions that this paper sought to answer. Our analysis demonstrates that, for M = 2700 GeV and / = 3.3, 
no special tuning or choice of parameters is required at the M\ scale to achieve B — L symmetry breakdown. 
Further analysis-for other values of M and /-shows that the percentage of points that break B — L is, in 
general, independent of M. This is because B — L breaking is not a question of generating a specific scale. 
Rather, it involves having soft masses aligned in such a way as to allow the S'-terms to drive the third family 
right-handed sneutrino tachyonic, see Eq. (58). On the other hand, B — L breaking is dependent on the 
choice of /. In the limit / —> 1, Sbl,r —> 0 and B — L breaking becomes impossible-whereas increasing / 
will allow B — L breaking to occur. A second feature of Table III is that, for M = 2700 GeV and / = 3.3, 
B — L breakdown where the Zr mass exceeds the experimental lower bound, although less prevalent, is 
still rather robust. In contrast to arbitrary B — L breaking, however, since Mz R — V2\m,y \ it follows that 
passing the Zr mass bound check is sensitive to the choice of M-with the survival rate increasing with M. 

A third important conclusion drawn from Table III is that, for the main scan, a large percentage of the 
initial points that have B — L breaking consistent with the Zr mass lower bound, also lead to the radiative 
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breakdown of electroweak symmetry. Note that the Z mass can always be adjusted-albeit by low energy 
tuning-to its experimental value of 91.2 GeV. Again, further analysis-for other values of M and /-shows 
that electroweak breaking, like B — L breaking, is also roughly independent of M. However, unlike B — L, 
small / favors electroweak symmetry breaking. As is well known from the literature [73-75], electroweak 
breaking occurs for universal boundary conditions-that is, for / = 1. On the other hand, as / increases, 
the randomly generated parameter m 2 Hu (Mi) can be considerably larger than the square of the initial stop 
masses. In this case, the RGE evolution to the SUSY scale may be insufficient to render m 2 H ^ (Msusy) 
tachyonic. However, since the initial soft masses arc randomly generated, the electroweak breaking survival 
rate will decrease with increasing /, but will not go to zero. 

Whether or not stop masses remain non-tachyonic at low scale depends on the randomly chosen values 
of several of the initial parameters. As can be seen from Table III, for the values of M and / chosen for 
the main scan, non-tachyonic stops arc very common. To remind the reader, the checks labeled spill arc 
repeats of earlier checks that are conducted after the final iteration of solving for the SUSY scale. Since 
this iterative process usually only affects the relevant checks logarithmically, the spill bins arc expected 
to have high survival rates. The survival rate for convergence of the iterative process of finding values of 
Mb~l and Msusy is almost 100%-since the soft masses have a logarithmic dependence on the scale. The 
survival rate for the SUSY particle mass bounds check is, for M = 2700 GeV and / = 3.3, comparable 
to that of the Zr mass bound. Further analysis shows that this rate is also controlled by the choice of M-a 
higher value for M resulting in a higher survival rate for this check. The Higgs mass survival rate for the 
main scan is, perhaps, surprisingly high-given that we were checking that the Higgs mass for a randomly 
generated point matches an experimentally measured value within an error of less than one percent. The 
reason this rate is so high is that the measured value of the Higgs mass is within the range expected for TeV 
scale supersymmetry breaking. 

Since the initial soft SUSY breaking parameter space is 24-dimensional, graphically displaying the 
subspaces associated with each survival check in Table III is very difficult. However, as can be seen from 
the RGEs and has been discussed in the text, much of the scaling behavior of the parameters is controlled by 
the two S'-terms, Sr and Sbl , defined in Eqs. (A24) and (A23). It follows that the results in Table III can be 
reasonably displayed in the two-dimensional Sbl(M\) - Sr{Mi) plane. We begin by presenting in Fig. 5 
the initial points in the Srl(M\) - Sr(M\) plane that satisfy, sequentially, the first two fundamental checks 
in Table III; that is, B — L breaking and the experimental Zr mass lower bound. Points that do not break 
B — L are shown in red, points that satisfy B — L breaking but not the Zr mass bound are in yellow, and 
points that break B — L symmetry and satisfy the Zr mass bound are shown in green. This plot reaffirms 
the conclusion drawn from Table III that B — L breaking consistent with present experiments is a robust 
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no B— L breaking 
M y bound violated 


My bound satisfied 


FIG. 5. Points from the main scan in the Sbl(M\) - Sn(M\) plane. Red indicates no B — L breaking, in the yellow 
region B — L is broken but the Zn mass is not above its 2.5 TeV lower bound, while green points have M Zft above 
this bound. The figure expresses the fact that, despite there being 24 parameters at the UV scale scanned in our work, 
B — L physics is essentially dependent on only two combinations of them-the two 5-terms. Note that the green 
points obscure a small density of the yellow and red points behind them. Similarly the yellow points obscure some 
red points. 


phenomena. Furthermore, it shows the strong dependence of B — L breaking and the Zr mass on the values 
of the .S'-tcrms. There is a line in the Sbl - Sr plane-between the yellow and red regions-below which 
B — L breaking is not possible. Note that this includes the origin, which corresponds to vanishing 5-terms 
and, hence, universal soft masses. This shows that at least a small splitting from sfermion universality is 
required for B — L breaking. Another line exists-between the green and yellow regions-below which Zr 
is always lighter than its lower bound. 

Proceeding sequentially, we present in Fig. 6 the initial points in the Srl{M{) - Sr(M{) plane that, 
in addition to breaking B — L with a Zr mass above the experimental bound, also break EW symmetry. 
The entire colored region encompasses the green points shown in Fig. 5. Those points that also break EW 
symmetry arc displayed in purple. This plot reaffirms the conclusion drawn from Table III that most of the 
points that break B — L with a Zr mass above the experimental bound, also break EW symmetry. Note that 
a small density of green points that do not break EW symmetry are obscured by the purple points. 

In Fig. 7, we reproduce Fig. 6 but now, in addition, sequentially indicate the points that are consistent 
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FIG. 6. A plot encompassing the green region in Fig 5. The green points in this plot correspond to those which 
appropriately break B — L symmetry, but which do not break electroweak symmetry. However, the purple points, 
in addition to breaking B — L symmetry with an appropriate Zr mass, also break EW symmetry. Note that a small 
density of green points that do not break EW symmetry are obscured by the purple points. 


with the remaining checks-that is, non-tachyonic stops/spill checks/convergence/all lower bounds on spar- 
ticles masses satisfied and, finally, that they reproduce the Fliggs mass within the experimental uncertainty. 
Points that appropriately break B — L symmetry but do not satisfy electroweak symmetry breaking are still 
shown in green. Points that, additionally, do break electroweak symmetry are again shown in purple. Such 
points that also have non-tachyonic stops, pass all spill checks and convergence, and which satisfy all lower 
bounds on sparticles masses, but do not match the known Fliggs mass, are now indicated in cyan. Finally, 
points that satisfy all checks, including the correct Fliggs mass, arc shown in black. These are the valid 
points. The density of black points indicates, as observed above, that there is a surprisingly high number 
of initial parameters that satisfy all present low energy experimental constraints. The distribution of black 
points can be explained from the fact that, while B — L breaking favors non-zero ,S'-tcrms. very large 5- 
terms can effect the RGE evolution of sfermion masses adversely. Since the effect of the 5-terms depends 
on the charge of the sfermion in question, some sfermions will become quite heavy while others light or 
tachyonic. Therefore, in general, the valid points in our scan are a compromise between large 5-terms, 
needed for a Zr mass above its lower bound, and small 5-terms needed to keep the sfermion RGEs under 
control. 
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FIG. 7. A plot of the valid points in our main scan. The green and purple points correspond to the green and purple 
points in Fig 6. The cyan points additionally have non-tachyonic stops, pass all spill checks and convergence, and 
satisfy all sparticle mass lower bounds. The black points are fully valid. That means that, in addition to satisfying 
all previous checks, they reproduce the correct Higgs mass within the stated tolerance. The distribution of points 
indicates that while B — L breaking prefers large S'-terms, sfermion mass constraints prefer them to be not too large. 
Again, the cyan and black points may obscure a low density of other points not satisfying their constraint. 


The most important property of the inital SUSY parameter space in determining low-energy phe¬ 
nomenology is the identity of the LSP. Recall that when /i-parity is violated, no restrictions exist on the 
identity of the LSP; for example, it can carry color or electric charge. Our main scan provides an excellent 
opportunity to examine the possible LSP’s and the probability of their occurence . To this end, a histogram 
of possible LSP’s is presented in Fig. 8-with the possible LSP’s indicated along the horizontal axis, and 
log 10 of the number of valid points with a given LSP on the vertical axis. The notation here is a bit con¬ 
densed, but is specified in more detail in Table IV. The notation is devised to highlight the phenomenology 
of the different LSP’s, specifically then - decays, which are also presented in Table IV. 

The most common LSP in our main scan is the lightest neutralino, X\- However, not all x? states are 
created equal. LHC production modes for the lightest neutralino depend significantly on the composition 
of the neutralino-a bino LSP cannot be directly produced at the LHC, but the other neutralino LSP’s can. 
This is the basis we use for the division of these states. The state designates a mostly rino or mostly 
blino neutralino, a mostly wino neutralino and \\ a mostly Higgsino neutralino. Here, the word mostly 
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LSP 


FIG. 8. A histogram of the LSP’s in the main scan showing the percentage of valid points with a given LSP. Sparticles 
which did not appear as LSP’s are omitted. The y-axis has a log scale. The dominant contribution comes from the 
lightest neutralino, as one might expect. The notation for the various states, as well as their most likely decay products, 
are given in Table IV. Note that we have combined left-handed first and second generation sneutrinos into one bin, 
and that each generation makes up about 50% of the LSP’s. The same is true for the first and second generation 
right-handed sleptons and sneutrinos. 

indicates the greatest contribution to that state. As an unrealistic example, if Xi is 34% wino, 33% bino 
and 33% Higgsino, it is still labeled . The chargino LSP’s are similarly separated into wino-like and 
higgsino like charginos, and the stops and sbottom divisions are as in our earlier papers, references [31, 32]. 
Note that this notation for the stops, t a d and i r , are only used to describe stop LSP’s. For non-LSP stops, 
we use the conventional notation t\ and G. 

To make Fig. 8 more readable, we have made an effort to combine bins that have similar characteristics. 
The first and second generation left-handed sneutrinos are combined into one bin, where about 50% of 
the LSP’s are first generation sneutrinos. The same holds true for the first and second generation right- 
handed sleptons, while the first generation right-handed sneutrino is always chosen to be lighter than the 
second generation right-handed sneutrino. This similarity between the first and second generation sleptons 
is expected, since their corresponding Yukawa couplings are not large enough to distinguish them through 
the RG evolution. For both sleptons and squarks, more LSP’s exist for the third generation-as expected 
from the effects of the third generation Yukawa couplings, which tend to decrease sfermion masses in RGE 
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Symbol 

Description 

Decay 

Y°~ 

a B 

Y°~ 

A w 

Xuc 

■y°_ 

X H 

A bino-like neutralino, mostly rino (Wr) or mostly blino (//')• 

Mostly wino neutralino. 

Mostly third generation right-handed neutrino. 

Mostly Higgsino neutralino. 

e±w*, vz, vh 

X* 

Xfj 

Mostly wino charginos. 

Mostly Higgsino charginos. 

uw ± , e ± z, e ± h 

9 

Gluino. 

ttv , tb£~ 

tad 

Left- and right-handed stop admixture. 

t+b 

t r 

Mostly right-handed stop (over 99%). 

tv, T + b 

9R 

Right-handed first and second generation squarks. 

£ + j, vj 

b L 

Mostly left-handed sbottom. 

bv 

t>R 

Mostly right-handed sbottom. 

bv, £~t 

^1,2 

Pl 3 

First and second generation left-handed sneutrinos. 

LSP’s are split evenly among these two generations. 

Third generation left-handed sneutrino. 

bb, W+W~, ZZ, 

tt, £' + t ~, hh, vv 

LRl,2 

First and second generation right-handed sneutrinos. 

vv 

tl 

Third generation left-handed stau. 

tb, W~h, 

ev, /MS, TIS 

cr, 9r 

First and second generation right-handed sleptons. 

LSP’s are split evenly between these two generations. 

eu, /ms 

tr 

Third generation right-handed stau. 

tb, ev, /j,v, tv 


TABLE IV. The notation used for the states in Fig. 8 and their probable decays. More decays are possible in certain 
situations depending on what is kinematically possible and the parameter space. Gluino decays are especially depen¬ 
dent on the NLSP, here assumed to be a neutralino. Here, the word “mostly” means it is the greatest contribution to 
the state. The symbol t represents any generation of charged leptons. The left-handed sneutrino decay into t' + t~ 
indicates a lepton flavor violating decay-that is, l' + and t~ do not have the same flavor. Note that j is a jet-indicating 
a light quark. 

evolution. 

The myriad of possible LSP’s leads to a rich collider phenomenology. This phenomenology is not 
the main focus of this paper, but it is worthwhile to briefly review it here. In models where /(-parity is 
parameterized by bilinear /(-parity, such as the B — L MSSM, SUSY particles are still pair produced and 
cascade decay to the LSP. At this point, the bilinear /(-parity violating terms allow the LSP to decay. While 
only a few studies have been done on the phenomenology of the minimal B — L MSSM [17, 18, 31, 32], 




35 


there have been several works on the phenomenology of explicit bilinear 77-parity violation, which has some 
similarities to this model. References to such papers arc mentioned below although see [76-79] for general 
discussions. Table IV provides some basic information on the most probable decay modes of each of the 
possible LSP’s. Note that i signifies a charged lepton of any generation and j a jet-implying a light quark. 
Interesting aspects of Table IV arc the following. 

LSP Phenomenology 

• Neutralinos: Only neutralinos with non-significant blino or rino components can be significantly 
produced at the LHC. Note that in addition to the usual possibilities, a mostly right-handed third 
generation neutrino is also a possible lightest neutralino component here, because of 77-parity vi¬ 
olation. This can be pair produced through the Zr resonance. Due to the Majorana nature of the 
neutralinos, they can lead to same-sign dilepton signals-a clear sign of lepton number violation. 
This is true whether they arc directly produced or occur at the end of a cascade decay. The genera¬ 
tion of l depends on the neutrino mass hierarchy, as discussed in [31, 32]. In the normal hierarchy, 
muons and taus are most likely, while in the inverted hierarchy all charged leptons are possible. 

• Gluino: Most of the LSP decay products mentioned in Table IV mimic well-known hypothetical 
states-for example, neutralinos decay like TeV scale right-handed neutrinos and squarks decay like 
leptoquarks. The same can not be said of the gluino, making it an interesting candidate for further 
study. However, its decays depend strongly on the identity of the next to lightest supersymmetric 
particle (NLSP). Also, bounds on the gluino arc the strongest because of its large production cross 
section. Therefore, when the gluino is the LSP, it is likely that it is the only LHC-accessible SUSY 
particle. As with the neutralinos, the gluino’s Majorana nature allows same-sign dilepton final states- 
indicating lepton number violation. 

• Squarks: All squark LSP’s act like leptoquarks in this model, meaning they arc pair-produced and 
decay into a lepton and a quark. Stop and sbottom LSP’s in this model were studied in [31, 32]. For 
both a down- and an up-type non-third generation squark LSP, there will be two highly degenerate 
LSP states-either a degenerate down and strange squark pair or a degenerate up and charm squark 
pair-as required by phenomenology. In the inverted hierarchy, these can decay into an electron and 
jet or a neutrino and a jet, making them tempting explanations for the recent CMS excess, see [80], 
in the eejj and evjj channels [81]. However, the branching ratios seem to be inconsistent with the 
cross section [82, 83]. See reference [84] for a study of stop in trilinear R-parity violation. 

• Left-handed sneutrinos: Left-handed sneutrinos decay like heavier neutral Higgs bosons, that is, 
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H° and ,4°, due to their /?-parity violating mixing with the Higgs sector. In general, decays into 
heavier Higgses arc also possible but, of course, this depends on kinematics. The final state t' + 1 
represent a lepton flavor violating final state, such as /x + e“. Sneutrinos LSP decays were studied in 
the case of explicit bilinear /^-parity violation, which has some similarities to the B — L MSSM, in 
reference [85, 86]. 

• Right-handed sneutrinos: These states decay into missing energy and, therefore, cannot be eas¬ 
ily distinguished from the /?-parity' conserving MSSM. However, since the sneutrino is spin 0, as 
opposed to spin half neutralinos, a detailed collider study might reveal some differences. It is also 
interesting to note that it may be possible to pair produce right-handed sneutrinos through a Zr 
resonance, although the cross section would probably be small. 

• Sleptons: Both left-handed and right-handed charged sleptons decay like charged Higgs bosons, 
with which the sleptons mix due to /f-parity' violation. The left-handed sleptons have more channels 
open to them because of their isospin charge. Each left-handed slepton comes in an SU (2) doublet 
with the associated left-handed sneutrino. Splitting of this doublet is mainly due to electroweak 
//-term contributions to the mass, which push the associated left-handed sneutrino to lighter mass 
values, making it the LSP. In the case of the left-handed stau, however, mixing effects through the 
Yukawa and tri-scalar couplings (see Appendix B 1) have the potential to make its mass lighter than 
the third-family left-handed sneutrino. Therefore, the left-handed stau is the only left-handed charged 
slepton capable of being the LSP. Slepton LSPs with explicit /(-parity violation were discussed in 
reference [87]. 

To get a sense of the non-LSP spectrum, we produce histograms of the masses of the sparticles from 
the main scan. In the following histograms, there will be quite a few pairs of fields that will be highly 
degenerate; these will be represented by only one curve. This includes SU{2)l sfermion partners, which arc 
only split by small electroweak terms. Lirst generation squarks arc also degenerate with second generation 
squarks with the same isospin, due to phenomenological constraints. A consequence of this is that all 
first and second generation left-handed squarks arc highly degenerate. In viewing these histograms, it is 
helpful to remember that aside from the usual RGE effects of the MSSM, there arc two additional effects 
involved. The first of these is the boundary conditions at the B — L scale, corresponding to the B — L and 
/;>/,< //-terms which arc given in Eq. (64). The second is the new RGE effects of the Sr and Srl terms. 
Although the signs of these terms are not fixed. Fig. 5 shows that Srl is typically positive while Sr is 
typically negative. This indicates that Sr will tend to increase (decrease) sfermion masses for sfermions 
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with a positive (negative) charge, while Sbl tends to increase (decrease) sfermion masses for sfermions 
with negative (positive) B — L. 

Figure 9 shows histograms of the squark masses. Because they come in SU (2) doublets and the first- 
and second-family squarks must be degenerate, all four of the first- and second-family left-handed squarks 
have nearly identical mass and the histograms coincide. The degeneracy of first- and second-family squarks 
is also evident in the right-handed squark masses. The first and second family right-handed down squarks 
are generally lighter than their up counterparts because of the effect of the C/( 1 ) 3/2 charge in the RGEs. 



mass(TeV) 


FIG. 9. Histograms of the squark masses from the main scan. The first- and second-family left-handed squarks are 
shown in the top-left panel. Because they come in SU (2) doublets, and the first- and second-family squarks must be 
degenerate, all four of these squarks have nearly identical mass and the histograms coincide. The first- and second- 
family right-handed squarks are shown in the top-right panel. The right-handed down squarks are generally lighter 
than their up counterparts because of the effect of the C/(1) 3 j?. charge in the RGEs. The third family squarks are shown 
in the bottom panel. 


Figure 10 shows histograms of the masses of the sneutrinos and sleptons. The third-family sleptons and 
left-handed sneutrinos tend to be the lighter because of the influence of the tau Yukawa couplings. The 
right-handed sneutrinos are labeled such that ub 1 is always lighter than Fr 2 . 
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FIG. 10. Histograms of the sneutrino and slepton masses in the the main scan. First- and second-family entries are in 
the top-left panel, along with the third family left-handed sneutrino. Staus are in the top-right panel with mass-ordered 
labeling. In the bottom panel, the first- and second-family right-handed sneutrinos are labeled such that D/ri is always 
lighter than i>R 2 - 

Figure 11 presents histograms of the CP-even component of the third generation right-handed sneutrino, 
the heavy Fliggses, the neutralinos, the charginos, and the gluino. The CP-even component of the third 
generation right-handed sneutrino is degenerate with Zr. It is always heavier than 2.5 TeV because we 
have imposed the collider bound on Zr. The neutralinos and charginos are labeled from lightest to heaviest 
as is canonical in SUSY models. The X 5 and are typically Higgsinos. 

We emphasize that all of the above histograms are calculated using our main scan; that is, for the choice 
of M = 2700 GeV and / = 3.3. We remind the reader that these values were chosen so as to maximize 
the number of valid accessible points. However, the mass scale of these histograms is heavily dependent on 
the choice of M. Smaller values for M will move the above distributions distinctly toward lighter sparticle 
masses. 

Plots of the physical particle spectra for four valid points are presented in Figs. 12 and 13. These four 
points are automatically selected from the pool of valid points from the main scan based on simple criteria. 
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FIG. 11. The CP-even component of the third-family right-handed sneutrino, heavy Higgses, neutralinos, charginos 
and the gluino in the valid points from our main scan. The CP-even component of the third generation right-handed 
sneutrino is degenerate with Zr. The and are typically Higgsinos. 


The first is the spectrum with an admixture stop LSP with the largest gap between stop LSP and the next 
lightest sparticle. The second is similar; now, however, with a right-handed sbottom LSP. The third and 
fourth are the valid points with the largest right-side-up and upside-down hierarchy respectively; that is, the 
largest splittings between the B — L and SUSY scales in the two possible hierarchies. 

Plots of the high-scale boundary values for four sample valid points from our main scan are presented in 
Figs. 14 and 15. While these look like Figs. 12 and 13, they do not correspond to physical masses but, rather, 
mass parameters at M\. These four valid points are automatically selected from the pool of valid points from 
the main scan based on simple criteria. The first two are those with the lightest and heaviest initial value 
of the third-family right-handed sneutrino mass. These show that it is not necessary to artificially choose a 
very light initial mass for the third-family right-handed sneutrino to effect the destabilizing of its potential 
and B — L symmetry breaking. Note when reading these plots that the lightest right-handed sneutrino is 
always, without loss of generality, defined to be the third-family. The next two plots show the valid points 
with the largest and smallest amount of splitting in the initial values of the scalar soft mass parameters. 
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FIG. 12. Two sample physical spectra with an admixture stop LSP and right-handed sbottom LSP. The B — L scale 
is represented by a black dot-dash-dot line. The SUSY scale is represented by a black dashed line. The electroweak 
scale is represented by a solid black line. The label up is actually labeling the nearly degenerate up and cp masses. 
The labels ur, dp and dp are similarly labeling the nearly degenerate first- and second- family masses. 




FIG. 13. Two sample physical spectra with a right-side-up hierarchy and upside-down hierarchy. The B — L scale 
is represented by a black dot-dash-dot line. The SUSY scale is represented by a black dashed line. The electroweak 
scale is represented by a solid black line. The label up is actually labeling the nearly degenerate up and Cp masses. 
The labels ur, dp and dpi are similarly labeling the nearly degenerate first- and second- family masses. 
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The amount of splitting is defined as the standard deviation of the initial values of the 20 scalar soft mass 
parameters. 




FIG. 14. Example high-scale boundary conditions for the two valid points with the lightest and heaviest initial value 
of the third-family right-handed sneutrino soft mass. The label Q i is actually labeling the nearly degenerate Q \ and 
Q 2 soft masses. The labels u c and d c are similarly labeling the nearly degenerate first and second family masses. 
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FIG. 15. Example high-scale boundary conditions for the two valid points with the largest and smallest amount of 
splitting. The label () \ is actually labeling the nearly degenerate () \ and Q 2 soft masses. The labels u c and d c are 
similarly labeling the nearly degenerate first and second family masses. 
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VII. FINE-TUNING 


Fine-tuning in supersymmetric models arises from Eq. (15), 

rn 2 H i tan 2 /3 — m 2 H 


-Ml = 


1 — tan 2 B 


l l ■ 


(77) 


In both the MSSM and the minimal B — L extension of the MSSM, the soft masses ml, and mi, receive 
contributions from other soft masses. Most important are the contributions from stop and gluino soft masses 
that appeal - in the RGEs for m 2 Hu ~see Eq. (A42). They must be TeV-scale to satisfy sparticle mass lower 
bounds and the measured Fliggs mass. These large TeV-scale contributions must be almost exactly cancelled 
to yield the relatively small value of M| on the left side of Eq. (77). This cancellation can come either 
from other soft masses or the fi 2 term. The delicate cancellation between the parameters on the right side to 
yield the smaller term on the left side is “fine-tuning”. The necessity of such fine-tuning in supersymmetric 
models has been referred to as the “little hierarchy problem”. 

Here we explain the little hierarchy problem within the context of the B — L MSSM, using a rough 
analytic argument along the lines of that presented in [88]. Although we discuss it using the language and 
notation of the minimal B — L extension of the MSSM, the same argument holds in the MSSM. The largest 
contributions to m 2 Hu come through its RGE in the B — L MSSM scaling regime, Eq. (A42). Focusing 
on just the stop and gluino soft mass contributions, we can write a solution to this equation to first-order in 
1ii(M[/Msusy)- Such a solution is quantitatively innacurate because it neglects higher powers of the large 
logarithm ln(Mi/MsusY)- Be that as it may, it can still provide insight into how various scales enter the 
problem. The solution is 


= -i r * 2( ”4+>4) + 


(78) 


where the ellipsis represents neglected higher order terms and terms due to other contributions in Eq. (A42). 
Additionally, there are corrections due to the boundary condition Eq. (64). The m 2 ^ and m 2 themselves 
receive large contributions through their RGEs, Eqs. (A44) and (A46). Focusing on the contributions from 
the gluino mass yields 


m H u = 


167T 2 


17 2 


"4 + -4 + In (VV) ) In (A4) + .... . (79) 


As discussed in Section IV A, the stops and gluino have relatively high mass bounds from LHC searches. 
Additionally, as discussed in Appendix B 2, satisfying the observed value of the Higgs mass tends to require 
heavy stops. This means that the stop and gluino soft mass contributions in Eq. (79) must be relatively 
large and give large contributions to the right-hand side of Eq. (77). For example, if Y t = 0.9, g 2 = 1, 
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Af> = m,q 3 = mjr = 1 TeV, M\ = 10 15 GeV, and Msusy = 1 TeV, these contributions are approximately 
equal to —(2 TeV) 2 . This must be almost exactly cancelled to yield the relatively small value of M§ = 
(91.2 GeV) 2 on the left-hand side of Eq. (77). The cancellation usually comes from the /r 2 , but can also 
come from the terms in the ellipsis or m 2 / ; . 

As stated above, in the minimal B — L extension of the MSSM there arc additional contributions coming 
from the boundary condition on the Higgs soft masses at the B — L scale, Eq. (64). Since tan j3 > 1, the 
most important arc the contributions to the H u soft mass. These arc proportional to Mz R . Since there is 
a lower bound of 2.5 TeV on Mz R , it is reasonable to suspect that these contributions to the H u soft mass 
necessitate more delicate cancellation-thus worsening the little hierarchy problem. Rewritten in terms of 
Mz r , the associated boundary condition is 

1 Mf„ . (80) 

z 9r t- 9 B l 


The gauge couplings here, and in the remainder of this Section, arc evaluated at Mb~l unless otherwise 
specified. Before concluding that this exacerbates the fine-tuning problem, we should replace the physical 
mass, Mz r , with more fundamental parameters of the theory, such as the soft masses evaluated at the 
intermediate scale. All of the scalar soft masses share in the generation of My n through the 5-terms. 
Substituting Eq. (61) into Eq. (17) allows us to write the 5-term contribution to Mz R . It is given by 

1 9r - 9r{ m i) 


My = - 


7 sfk(M. i) 

1 9bl ~ 9bl(M\) 


Sr(M{) 

Sb-l(Mi) + 


(81) 


^ 9bl(Mi) 

Substituting this into equation Eq. (80) yields 5-term contributions to the H u soft mass. In addition, 
the 5-terms also influence the running of the H u soft mass through the RGEs. Including both of these 


contributions, the value of rn 2 ,, at the SUSY scale is 

-tlu 


m H u { M SUSY) = 


9r 


1 9r 9r{Mi) 1 g 2 BL g BL (Mi) 

Sr(M,) “ 06 9 b L (Mi) Sb ~ l(M ' ] 


9r + 9bl V 14 9 2 r{M\) 

1 9 2 r-9 2 r(M i) 5fl(Ml) + 


(82) 


14 4(Mj) 

Consider a sample valid point with Sr(M\) < 0 and Sr-l(Mi) = —2Sr(M\). This case fits within the 
valid black points in Fig. 7. This arises physically if all scalar soft masses arc universal with the exception 
that the the first- and second-family right-handed sneutrino soft masses-which are heavier. In this case, the 
5-term contributions to rn 2 u can be written as 


m2 H u ( M SUSY) = 


9r 


1 9r 9r{Mi) 1 9rl 9bl(M i) \ 
9r + ^9bl V 44 9r{M\) 8 g BL (M\) ) 


1 9r ~ 9r(Mi) 
14 4(Mj) 


Sr(Mi) + 


(83) 
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Let us choose, for example, Msusy = 1 TcV and Mb~l = 2.5 TeV. Then the dimensionless coefficient of 
Sr{M i) turns out to be —0.022. Since this value is considerably smaller than unity, it follows that the //,, 
soft mass is not-in fact-very sensitive to the fundamental parameters that set the Zr mass. This remains 
true for all values of Msusy and Mb~l associated with valid points. Therefore, there is not a significant 
amount of new fine-tuning introduced in this way. 

Fine-tuning in supersymmetric models has historically [89-92] been quantified using the Barbieri- 
Giudice (BG) sensitivity, introduced in [93] and [94]. This quantifies the sensitivity of some observable 
quantity to changes in any of the fundamental parameters of a theory. The delicate cancellation between 
TeV-scale supersymmetry parameters in Eq. (77) results in the electroweak scale, Mz, having a large BG 
sensitivity. The BG sensitivity of the electroweak scale is defined as 


ai gMf 

ai M\ Qa, ’ 


(84) 


where a, is any of the fundamental parameters of the theory. This says that a fractional change in a, would 
produce a fractional change in M| that is F ai times larger. The overall degree of fine-tuning is usually taken 
to be the largest of all the F ai ’s ; that is, 


F = max(F ai ) . 


(85) 


The BG sensitivity F will be used to quantify the fine-tuning required in the B — L MSSM. 

It is worth mentioning that some authors have pointed out drawbacks to the BG sensitivy and suggested 
other quantifications of fine-tuning. For example, as discussed in [91], the BG sensitivity and the overall 
degree of fine-tuning depend on how the fundamental parameters of the theory a, are chosen. Furthermore, 
one could reasonably use the BG sensitivity of Mz, rather than M§, as the indicator of fine-tuning. This 
would result in fine-tuning that is smaller by a factor of two. Such ambiguities in the way fine-tuning is 
calculated from the BG sensitivity suggest that it is not a precise way to quantify fine-tuning. A separate 
paper, [95], points out that the relationship between the proton mass and the strong coupling constant at 
a high scale exhibits high BG sensitivity, but is not actually finely tuned. They propose a more precise 
quantification of fine-tuning and show that the BG sensitivity actually overestimates the tine-tuning in some 
sample points in the MSSM. 

Despite the possible shortcomings, the BG sensitivity remains the most widely used tool for making 
rough quantitative analyses of fine-tuning in supersymmetric models. We, therefore, proceed using the BG 
sensitivity to quantify tine-tuning in the B — L MSSM. For each of the valid points, we compute F. We 
allow a, to span all of the soft mass parameters of the theory, as well as //. In the case of scalar soft masses, 
we take a, to be the mass squared, while in the case of gaugino soft masses and / t we take a, to be the mass 
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to the first power. This choice corresponds to how these parameters appear in the Lagrangian. We then 
create a histogram of F for all of the valid points in our main scan. This data is shown as the blue line in 
Fig. 16. Note that the fine-tuning required by the highest percentage of valid points is F ~ 5000. Be that 
as it may, a reasonable number of valid points need significantly less fine-tuning-with about 2% requiring 
F < 2000. It is interesting to compare the amount of fine-tuning in the minimal B — L extension of the 
MSSM model to the amount of fine-tuning required in an identical statistical scan of the /(’-parity invariant 
MSSM using M = 2700 GeV and / = 3.3. Due to the aforementioned ambiguities in how fine-tuning 
is quantified, it is critical that the fine-tuning be calculated the same way when two different models are 
being compared. Therefore, we use our own code, slightly modified, to produce a similar plot for the 77- 
parity conserving MSSM. The results are shown as the green line in Fig. 16. Comparison of the blue and 
green lines in the figure show that the B — L MSSM valid points tend to be slightly less finely tuned than 
valid points in the 77-parity conserving MSSM. The difference is large enough to be apparent in the figure. 
Flowever, due to the unresolved questions about how to properly quantify fine-tuning, we do not regard this 
difference between the B — L MSSM and the MSSM to be significant. 



FIG. 16. The blue line in the histogram shows the amount of fine-tuning required for valid points in the main scan of 
the B — L MSSM. Similarly, the green line specifies the amount of fine-tuning necessary for the valid points of the 
77-parity conserving MSSM-computed using the same statistical procedure as for the B — L MSSM with M = 2700 
GeV and / = 3.3. The B — L MSSM shows slightly less fine-tuning, on average, than the MSSM. 


With the fine-tuning of each randomly generated point in the B — L MSSM now quantified, we are 
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equipped to produce results for just the most natural points-that is, those requiring minimal fine-tuning. 
Figure 17 shows a histogram of the LSP’s for those points with F < 1000, corresponding to the least 
fine-tuned ~ 0.1% of points, from a larger scan of four hundred million points. We refer to these points as 
“natural” valid points. There arc three notable differences between Fig. 17 and Fig. 8. First, stop LSP’s 
are more common. This includes both admixture and mostly right-handed stop LSP’s. Stop LSP’s arc 
more common because heavy stops tend to cause fine-tuning, so low fine-tuning favors lighter stops and 
stop LSP’s. Second, sbottom LSP’s arc more common. This is due to the fact that first, both the stop and 
sbottom masses depend on the soft mass mjy., and second, because the right-handed stop and sbottom soft 
masses have similar terms in their RGE’s. These two facts imply that favoring light stops tends to favor 
light sbottoms as well. Third, Fig. 17 does not have the gluino LSP’s shown in Fig. 8, and it does have some 
Hr LSP’s not found in Fig. 8. However, the disappearance and appearance of these states in the F < 1000 
histogram is not statistically significant and, hence, these states can be ignored. The prevalence of stop and 
sbottom LSP’s is the only significant difference between the natural valid points and the valid points. Note 
that the physical implications of having a stop or sbottom LSP in the B — L MSSM have been studied in 
[31, 32]. 


VIII. CONCLUSION 

In this paper, we presented a novel approach for relating UV physics to TeV scale physics and applied 
this analysis to the minimal SUSY B — L model. This approach hypothesizes that all SUSY breaking pa¬ 
rameters are about an order of magnitude away from a characteristic SUSY breaking mass scale. Practically, 
this translates into conducting an analysis where all relevant soft SUSY mass parameters arc independently 
scanned over the same range at the UV scale, and then RG evolved to the TeV scale. This program lends 
itself especially well to the string realization of the minimal B — L MSSM model. However, our results 
are relevant for any high scale soft SUSY breaking minimal SUSY B — L model with gauge coupling 
unification. 

A central result of this work is the general region of initial parameter space that leads to radiative B — L 
symmetry breaking. While this depends on multiple parameters of the theory, it can be expressed in terms 
of the two S-parameters and is presented in this context in Fig. 5. A subsequent figure. Fig 7, shows how 
additional constraints, such as electroweak symmetry breaking and lower bounds on new sparticle masses, 
depend on the S-parameters. These two plots indicate that a significant amount of the initial parameter 
space leads to experimentally viable results. They arc followed by various spectrum graphs which show 
that acceptable spectra are relatively general and do not depend on a specific hierarchy of initial masses. 
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LSP 


FIG. 17. A histogram of the LSP’s for the “natural” valid points with F < 1000. Sparticles which did not appear 
as LSP’s are omitted. The y-axis has a log scale. The notation for the various states, as well as their most likely 
decay products, are given in Table IV. Note that the natural valid points favor stop and sbottom LSP’s more than the 
valid points presented in Fig. 8. Note that we have combined left-handed first and second generation sneutrinos into 
one bin and each generation makes up about 50% of the LSP’s. The same is true for the first and second generation 
right-handed sleptons and sneutrinos. 

The phenomenology of a given point at the LHC strongly depends on the identity of the LSP. Therefore, 
another central result of this paper is the calculation of the probability that a given SUSY particle can 
be the LSP. This was addressed in Fig. 8. As might be expected, a mostly bino neutralino is the most 
likely candidate. However, since binos cannot be directly produced at the LHC, signals associated with 
bino LSPs also depend on the rest of the SUSY spectrum. Therefore, an interesting future direction might 
be to investigate the phenomenology of mostly wino or Higgsino neutralinos. Mostly wino or Higgsino 
neutralinos can be directly produced at the LHC, independently of the rest of the SUSY spectrum, and 
have relatively large cross sections for colorless particles. The signals associated with different LSPs arc 
summarized in Table IV. 

Finally, the fine-tuning associated with this statistical scan was investigated. While it is not drastically 
different than the fine-tuning in the MSSM with a similar UV completion, one might think that the new mass 
scale associated with B — L breaking could introduce new contributions to fine-tuning. We showed that it 
does not. In fact, a given point in this model is typically less fine-tuned than a similar point in the MSSM. 
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In addition, we explored possible LSPs for points with fine-tuning better than one paid per thousand-in a 
way analogous to Fig. 8 . We found that stops and sbottoms become much more likely LSP candidates, as 
one might expect-see Fig. 17. The signals of stop and sbottom LSPs were discussed previously in [31, 32]. 
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Appendix A: Renormalization Group Equations 


This Appendix lists the RGEs used in this study. Most RGEs arc derived with the help of reference [42], 
unless otherwise stated. 

The RGEs for gauge couplings were presented in Section III, but are repeated here for completeness. 
The RGE for a general gauge coupling is 


d_ _1 _ _K_ 

dt a 2i r 


(Al) 


where t is the logarithm of the renormalization scale and the index a runs over the different gauge factors. 
The slope factors are different in each of the different scaling regimes: 


• Intermediate regime: 63 = 10, b 2 = 14, b R = 14, bs~L = 19. 

• B — L MSSM: b 3 = -3, b 2 = 1, b 3R = 7, b B _ L = 6 . 


• MSSM: 63 = -3, b 2 = 1, h = f. 

• Non-SUSY B - L: b 3 = -7, b 2 = -f, b 3R = f§, b B . L = f. 

• SM: 63 = —7, b 2 = -f, b 1 = % 


The gaugino soft mass RGE is 


Jt M ■ = 


baCX-aMa 
27T ’ 


(A2) 


where the b a arc the same slope factors given in Eqs. (37 - 40). It is helpful to observe that the gaugino mass 
renormalization group equation admits a rather compact analytic solution: 


M m _ M a {Mu) 
-M a (t) — 0^a\t 

OL\J 


(A3) 
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for all gaugino masses associated with SO( 10) and 


M i(t) = -ttt -r-ai(t), 

ai{M B -L) 


for the bino. 


There are three significant Yukawa couplings for RGE analysis: y t , yb and y T . In the SM scaling regime 
their RGEs can be found in [96], for example, and are given by 

j t y t = Y^2 yt (l^ y t ~ + 3 (% 2 + Vb) + y r~ 8y i ~ \al ~ (A5) 

j t Vb = ~ y ^ + 3 ( y * + y2 ) + £ ~ 8 si ~ \al - \v \) (A6) 

d _ 1 /^ 3 .2 , 0/2 , _. 2 \ , 2 9 2 9 2 A /A ox 


Jt VT = 16^ \2 y T + 3(% +y ^ + y r- 4 d2 - 4 91 ) ■ 

In the U (1) extended SM regime of the upside-down case, the Yukawa coupling RGEs arc 

j t y t = ^6^2 y* Q(% 2 - Vb) + 3 (Vt + Vb) + yl 
-8gl ~ 1 5a - 4 ^ - * <?Il) 

dt yb = I67r 2 ^ b (y2^ b ~ ^ y ^ 

o 2 9 2 3 2 1 2 \ 

_ %3 _ 4^2 _ 4 _ 4RBL ) 

= (§!/? +!!(!/? + * + !/? 

9 2 3 2 9 2 \ 

— ^92 ~ 4 9r ~ ~^9bl ) • 


(A10) 


The boundary condition at the B — L scale is trivial. At the SUSY scale, however, the boundary condition 
is nontrivial: 


yt(M, susy) = ^(Msusy) sin (3 
yb , r { M susy) = Y ^ T { M , susy) cos/T 


(All) 


The Yukawa couplings above the SUSY scale will be denoted by Y instead of y. This condition applies 
both in the upside-down case and in the right-side-up case. In the MSSM scaling regime of the right-side-up 
case the RGEs are 


S^ = ^( 3r U«f-3 9 ?-frf). 


(A12) 


(A13) 
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In the B — L MSSM scaling regime the RGEs arc 

Y = Y y < ( 6F < 2 + y2 - f 9 * - 392 -- 4") <A14) 

j t Yb = Ife 5 n ( 6y2 + y2 + y2 “ ¥®> ” %2 " ii 9 ® 1 _ 9 h) <A15) 

Y = ( 3y2 + 4y2 - - 5 9 “ - 9 l») ' <A16) 

The fact that these RGEs arc non-linear means that the analytic solutions arc much more cumbersome if 
they can be found at all. We use numerical integration techniques instead, yielding numerical values for 
the Yukawa couplings at any scale up to the intermediate scale, M\. These solutions will be subsequently 
used in the running of the soft tri-scalar couplings and some of the scalar soft masses because the RGEs of 
those parameters depend on the Yukawa couplings. The Yukawa couplings do not need to be evolved above 
the intermediate scale since the couplings that depend on them will not be evolved above the intermediate 
scale. 

Tri-linear couplings arc generated at the intermediate scale and evolved to the SUSY scale. Their RGEs 
in the B — L MSSM scaling regime arc 

d _ 1 ( a 2 , v 2 16 2 q 2 1 2 2 \ 

Jt at ~ 16^°* V 8 + Yb - Y 93 ~ 392 ~ 6 9bl ~ 9 Y) 

+ 10^2 ^1 ® a uYt + 2 Ybdb + —g 2 M$ 

+6.c/| M2 + -g^ L Ms-L + ‘^gjRMjii'j (A17) 

d 1 / 0-5^2 1 -\r2 , v 2 16 2 q 2 1 2 2 | 

Jt ab - wY ab + K + Yt ~ Y 93 ~ 392 ~ q 9bl ~ 9 Y) 

+ 10^2 Y (^-^ a b^b + 2 a T Y T + 2Yta u + —g 2 M% 

+6g|M2 + -g^ L Ms-L + ‘IgjrtMiR^ (A18) 

Jt ar = YY ar { 3Yb + 6Y " ~ 392 ~ \ 9bl ~ 9 l R ) 

+ W^Y t ^6 abYb + 6a T Y T + 6g^M2 + 2>g 2 ^^MB-L + (A 19) 


In the right-side-up case, the B — L scale is above the SUSY scale so these parameters will also be run 
through the MSSM scaling regime from the B — L scale to the SUSY scale. The RGEs in the MSSM 
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scaling regime arc 


d 


dt H 


d 


dt 


ab 



+ ^ ^10a u lt + 2Tf/U + — 53M3 + 6ry|AT, + 

+ 10^2 ^ ^10 a 6^6 + 2a T Y T + 2Yta u + —g^M^ 
+652 M2 + — gfMi'j 
^2 a r (W 6 2 + 6F t 2 - 3^1 - 
+ ^^2 (^ a b^b + ha T Y T + 6<7! M2 + —g 2 M\\ . 


(A20) 


(A21) 


(A22) 


These equations are also do not yield tractable analytic solutions, of course. 

Scalar soft mass squared parameters are also inputted at the intermediate scale and evolved down to 
the SUSY scale. In the case of the right-side-up hierarchy, this will involve running through the B — L 
scale and the brief MSSM scaling regime. The boundary condition at the B — L scale is nontrivial because 
U-term interactions between the third-family right-handed sneutrino and the other scalars give rise to a new 
contribution to the soft masses when the third-family right-handed sneutrino acquires a VEV, Eq. (64). As 
discussed in Section IV B, we take the soft masses to be flavor diagonal in order to satisfy flavor constraints. 

Before writing the scalar soft mass RGEs, it is useful to define the S'-terms, 

Sb-l = Tr (2mj~ - m| c — m| c — 2m| + m~ c + m~ c ) (A23) 

Sr = m 2 Hu - m 2 Hd + Tr + |m 2 ~ c - ^m| c + (A24) 

Sy = m 2 Hu - m 2 Hd + Tr (mj - 2m? c + m 2 ~ c + m\ - mf c ) , (A25) 


where the traces arc over generational indices. It can be shown, using the scalar soft mass RGEs, that the 
5-terms obey the RGEs: 


d q _ b a a a S a 

Jt a = 2-7T 


(A26) 


which admit the simple analytic solution 


for Sr or Sb-l and 




Sr(t)= s m^ Sr(MsmY> ’ 


(A27) 


(A28) 



52 


It is perhaps useful to separate the scalar mass RGEs into those that arc analytically tractable and those 
that are not. In the B — L MSSM scaling regime, the first- and second-family and sneutrino soft mass RGEs, 
analytically solvable, arc 


1671-2 1,2 = - y33 M 3 - 6<£ M 2 - ^ 9bl m b-l + \9blSb-l (A29) 

2 = -y 5 |Af| - g BL M B-L - 2 g 2 R M 2 R 

— ~^9blSb~l — QrSr (A30) 

= ~ Y ff 3 M 3 “ \9bl m b-l ~ 2 (Jn M n 

~~^9bl^b-l + 9rSr. (A31) 

16tE— m| i2 = — 6 r /2 — 3g BL M B _ L — -g^ L S B -L (A32) 

167T 2 —m|c 2 3 = —3 g B L^B-L ~ %g R M R + ~^9blSb-l ~ 9rSr. (A33) 

167T 2 —m|c 2 = —3 9b L M\_ l — 2g R M R + —g BB S B —L + 9 rSr- (A34) 


In the MSSM scaling regime, which is only relevant to the case of the right-side-up hierarchy, the RGEs arc 


ir 2 d 2 

l07T —m^ 

dt Q 1.2 

-ic 2 ^ 2 

l07T —m.-.c 

dt, u 1.2 

ic 2 ^ 2 

107T —mj r 

dt d 


1,2 


d 


-\c 2 u 2 

l07T —mf 

dt G,2 

1C 2 ^ 2 

107T —m,-.c 

dt, A,2 

1C 2 ^ 2 

lo7r —rrizc 
dt, A ,2 


-f 9?Mf - 6„lM| - + IjfSv 

-§, - gsfMf - 
-f 9i«l - A 9 2 Ml 2 + A//PV 
- A.»r - A 9 ?s k 

5 5 

0 

-|y 2 5 2 M 2 + |yp 2 5 r 

5 5 


(A35) 

(A36) 

(A37) 

(A38) 

(A39) 

(A40) 


The right-handed sneutrinos masses do not run in this regime because they arc not charged under the MSSM 
gauge group. In the upside-down case the the right-handed sneutrinos are present in the brief scaling regime 
between Msusy and M b ~l ■ Their soft mass RGEs arc 


1C 2 u 2 
l07T —m.-,c 

dt A,2,3 


° 2 / 2 , 2 , 2 \ 
T 9BL\ m vl + 


(A41) 


For the third family sfermions (excluding the sneutrinos) and for the MSSM Higgs, all of which arc not 



53 


analytically solvable, the RGEs In the B — L MSSM scaling regime are 


= 6Y t 2 ( m H u + m | 3 + m %) + 6a ? 

—6 p g|iVf| — 2 gjiMft + (JrSr 

—rnjjd = 6 Y${mJj d + m. 2 ^ + ra~ c ) + 2 Y 2 {m 2 Hd + m| g + rrif c) + 6a 2 + 2a 2 
—6g|Af| — 2 g^Mft — g\Sn 

= + m | 3 + m D + 2 Y b( m H d + + m b c ) + 2 «l + 2 «b 

-Y93 M 3 - %2 M 2 - \d 2 BL M BL + \gBL S B-L 

l67r2 J t m h = 2y ^ m ^d + m ! 3 + m f c ) + 2a r 

—6^|M| — 3 gBL^BL — '~^9 bl^b-l 

l ^ 2< lit m i c = 41 f( m k + m | 3 + m t=) + 4a ? 

~^~93 M 3 - 2 9bl m bl - 2 9r m r ~ 9bl s b-l ~ 9rSr 
= 4 Y b( m H d + m | 3 + m fc) + 4afe 

~^93 M 3 ~ g 9bl m bl ~ 2 9r m r ~ ^9bl s b-l + 

167r 2 -^m| c = 4Y 2 (m 2 Hd + m| g + m~ c ) + 4a 2 

~^9bl^bl ~ 2 9rMr + ^ 9bl^b-l + 9rSr- 


(A42) 

(A43) 

(A44) 

(A45) 

(A46) 

(A47) 

(A48) 


In the MSSM scaling regime they are 


1P 2 ^ 2 
lb7T —mu 

dt Hu 


in 2 ^ 2 

16 ' s m «- 


in 2^ 2 

167T — m.= 

dt Qz 


6 Y t 2 {m 2 Hu + m| 3 + m|) + 6 a f 2 

-6<7 2 2 M 2 2 - ^ 2 M 2 + 

5 5 

6 Y d{ m H d + m ^ 3 + m fc) + 2 Y 2 (m 2 Hd + m| 3 + mfc) + 6a 2 + 2a 2 

-6ff 2 2 M 2 - p 5 2 Af 2 + | 5 ?Sy 
5 5 

2 Y t{ m H u + m^ 3 + m~ c ) + 2 Y 2 (m 2 Hd + m| 3 + mg c ) + 2 a 2 + 2a 2 

-Y93 M 3 - 6 - 92 M 2 - Y^9\ M l +^gfS Y 


(A49) 


(A50) 


(A51) 
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= 2 Y r( m H d + m l 3 + m f «=) + 2a ? 

-6.g 2 2 M| - y5 ?M 2 - plS Y (A52) 

= 41^ 2 (jn 2 Hu + mj g + m fc ) + 4 a 2 t 

(A53) 

16 - 2 |-| c = 4y fe 2 (m# d + m| s + m? c ) + 4ag 

32 on 8 o , 2 o 

- y5 2 M| - - <? 2 M 2 - -gfSy (A54) 

167r 2 y m 2 c = 4y r 2 (m^ d + m| g + m 2 c ) + 4a 2 

-.B 5 2 M 2 + ^2 5y (A55) 

The soft mass parameters are used in the calculation of the physical sparticle masses, discussed in the next 
appendix. 

Appendix B: Physical Masses 

In this Appendix, we discuss how the physical masses of the sparticles and the Higgs are determined 
from the running parameters. 

1. Sparticle Masses 

Because the first- and second-family Yukawa and tri-scalar couplings are negligible, mixing among the 
first- and second-family sfermions and the sneutrinos is negligible, greatly simplifying the relationship be¬ 
tween physical masses and soft masses. However, there are electroweak 14-term contributions associated 
with the electroweak scale. Although these are numerically small, they have the effect of splitting the 
masses of the otherwise degenerate SU (2 )l doublets, which has implications for the lightest supersymmet¬ 
ric particle (see Section VI): 

= Mf (T 3 — Q sin 2 Ow) cos 2/3, (Bl) 

where 0\y is the weak mixing angle (sin 2 0\y ~ 0.23) and T : > and Q are the left-handed isospin and electric 
charge of the scalar (/>. Here we lay out the physical masses with the electroweak 14-term contributions, 
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along with the notation for the physical masses. 


ma L 

= mu 

+ 

A Qi’ 

m UR 

= 

rriu c 

+ A « c ) 

m ~c L 

= nric 

+ 

A Q 2 ’ 

m c R 

= 

rricc 

+ A 3 c, 

m d L 

= m d 

+ 

A <3i ’ 

m d R 

= 


+ A cZ> 

™S L 

= ms 

+ 

A (?2’ 

m~ SR 

= 

msc 

+ Age, 

™*L1 : 

= m,o x 

+ 

A Li’ 

m t>Rl 

= 

m D c 

+ A i>f: 

m VL2 ' 

= mo 2 

+ 

A l 2 ’ 

mo R 2 

= 

mpe 

+ A ^2' 

TOp L g : 

= mo 3 

+ 

A Z 3 > 

m 

V R 

= m Z r , 

m-e L 

= ms 

+ 

A Zi’ 

m~e R 

= 

mge 

+ Age, 

mfi L 

= mp, 

+ 

a l 2 ’ 

m »R 

= 

mp c 

+ A^c. 


The third-family right-handed sneutrino physical state (referred to as ifp) mass is different because it ac¬ 
quires mass through the B — L symmetry breaking mechanism and is degenerate with the Z/ ; mass. 

The Yukawa and tri-scalar couplings associated with third-family squarks and charged sleptons con¬ 
tribute non-negligible mixing terms among these scalars. These effects are captured in the stop, sbottom, 
and stau mixing matrices. Here we use the conventional notation a t ^,r = Yt,b, T ^t,b,T- The stop mixing 
matrix in the basis (t. t c *) is 


M\ = 


m\ + M? 
Qz 


A 


Qz 


M t ( A t - 


tan (3 


M t ( A t - tizp) m % + M t + 


(B3) 


The eigenstates of this matrix will be referred to as t \ and £2 with mass eigenvalues defined such that 
nip < mp 2 . The sbottom mixing matrix in the basis ( b , If*) is 

m| 3 + M 2 + A q 3 M b (A b - [i tan f3) 


Mi = 
0 


M b (. A b - n tan /3) m? c + M ft 2 + A~ b 


(B4) 


The eigenstates of this mass matrix will be referred to similarly to the stops. The stau mixing matrix in the 
basis (f, e c *) is 


M 2 = ( m l 3 + M t + A Z 3 M r ( A t ~ V tan /?) \ ^ (B5) 

y M r {A t — n tan (3) m^c + M 2 + Afc J 

The eigenstates of this matrix will be referred to similarly to the stops and sbottoms. All of the running 
parameters in these matrices are evaluated at the SUSY scale. 


We present these matrices in terms of the fermion masses Mtj>, T for simplicity. However, for numerical evaluation these fermion 
masses are replaced with the appropriate Higgs VEV times Yukawa coupling evaluated at the SUSY scale. 



56 


For any of these matrices, 


( L f X f ) 

V x f R f) 

the relevant mixing angle is given by 


(B6) 


-21 X f \ 

tan 29 ? = , (B7) 

L f ~ R f 

where the angle 9j may always be chosen to be between 0° and 90°. Defined this way, a mixing angle close 
to zero means the lighter mass eigenstate consists of mostly the left-handed gauge eigenstate and a mixing 
angle close to 90° means the lighter state is mostly right-handed. 

The chargino content is identical to that of the MSSM in the approximation of vanishing /i-parity vio¬ 
lation. This is a good approximation for calculating masses but the mixing with the charged leptons need to 
be take into account when calculating decays, see [31] for example. Continuing with the approximation of 
vanishing imparity violation, the results of [3] may be used. Those results, in our own notation, arc 


mt± — — (m$ + p 2 + 2 Myy — \J (Mf + [i 1 + 2Myy) 2 — 4(/xM2 — sin 2/3) 2 ^ 
m lf = 2 ( M 2 + I 1 ' 2 + 2M W + ^(Ml + » 2 + 2 - 4 GuM 2 - Myy sin 2p£\ . 

In the basis (v, Wr, B 1 , W°, Hy , H®), the neutralino mass matrix is 


l 

0 

c o r m z r 

s 6 r M Zr 

0 

0 

0 


-C 0 r M Zr 

Mr 

0 

0 

-cpsg w M z 

spsg w M z 


so r m z r 

0 

Mbl 

0 

0 

0 


0 

0 

0 

m 2 

C 0 Cg w M Z 

-SpCg w M Z 


0 

-cpsg w Mz 

0 

cpcg w M z 

0 

-F 

V 

0 

S 0 Sg w M Z 

0 

~SpCg w M Z 


0 


(B8) 

(B9) 


(BIO) 


where cq = cos 9 and sg = sin 9 etc. As with the charginos, we have assumed that mixing with the left- 
handed neutrinos, due to /7-parity violation is 0. This is good approximation for calculating masses and 
will be used here, but cannot be used when calculating decay rates. As discussed in Section IV some of the 
eigenstates of this matrix have masses associated with the B — L scale while others have masses associated 
with the SUSY scale. A conventional approach to this situation would be to perturbatively diagonalize the 
matrix in the limit A/.susy A> Mb~l for the right-side-up case or Mb~l A/susy for the upside-down 
case. However, these two scales may be comparable so the entire mass matrix must be diagonalized without 
the use of perturbative methods. This has the potential to introduce errors since it doesn’t account for the 
fact that some states should be integrated out at different scales. However, the errors will always be small 
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because the B — L and SUSY scales arc always of comparable size. We choose to evaluate all of the running 
parameters in this matrix at the SUSY scale. The error introduced by doing this should be smaller than the 
error introduced by associating the entire SUSY spectrum with a single scale, AfsusY- The mass eigenstates 
are referred to as Xi '' ■ Ac in a mass ordered basis with eigenvalues rri^o ■ ■ ■ m^o. 

The physical gluino mass, Mg is simply equal to the running gluino mass evaluated at the SUSY scale. 

Mg = M 3 ( Msusy ) • (B11) 


2 . Higgs Masses 


Supersymmetric models such as the MSSM and this B — L MSSM contain five Higgs particles. The 
most important for the present discussion is the lightest neutral SM-like Higgs, h°, which we refer to as 
“the Higgs” throughout this paper. This one is important because its mass is known and can be used to 
constrain some of the SUSY parameter space. The other four Higgses are the heavy Higgs, H°, the Higgs 
pseudoscalar, A °, and the charged Higgses, H ± . 

The Higgs mass is calculated using methods discussed in [46^18]. The physical Higgs mass is 

m h o = '/Xv, (B12) 


with the Higgs quartic coupling. A, evaluated at the scale of the physical Higgs mass. Above the SUSY 
scale, A comes from the D-terms and is thereby fixed. Below the SUSY scale, RGE effects will cause A to 
deviate from its supersymmetric value. These effects come mainly from one-loop graphs involving the top 
quark. They are contained in the RGE for A in the SM scaling regime. We employ results from [48]. Here 
we re-state the relevant equations in our own notation. The supersymmetric boundary condition on A is 

A(Msusy) = ^ ^9l + cos2 2 /3 + (B13) 


The parameter 5 A contains threshold corrections applied at the SUSY scale. Including only the dominant 
stop contributions from [48], 


167t 2 5A 


3Y t 4 


2 X * F 
m h m t2 



1 xf 

6 TO? TO? 
1 1 1 2 


G 



(B14) 


where we define X t = .4/ — // cot ft (note that this definition is different from that used in [48]) and 


F(x) 
G(x ) 


2x In x 
x 2 — 1 

12x 2 (1 — x 2 + (1 + x 2 ) lnx') 
(x 2 — l) 3 


(B15) 


(B16) 



58 


The RGE for A in the SM regime is 

—A = 4A(3yf + 3 y 2 + yf) — 9X(-gf + g 2 ) 

97 Q Q 

“4(3 yf + 3 y b + y T ) + ^qSi + <?25i + ~^gi + 12A 2 , (B17) 

and in the upside-down case between Msusy and Mb-l it is 

— A = 4A(3y 2 + 3 y 2 + yf) — 9A(-p|. + g 2 ) 

“4(3 yf + 3r/ft + yf) + -g 4 fi + -g 2 gji + ~^g\ + 12A 2 . (B18) 

Since this depends on the Yukawa couplings, which arc solved numerically, this must also be solved numer¬ 
ically. The dominant contributions come from the terms involving y t . These terms are present because both 
stops arc integrated out at Msusy- This has the potential to introduce errors because the stops generally do 
not have the same mass. The errors introduced by this arc minimized when the SUSY scale is chosen to be 
Msusy = J m t 1 m u ■ We find this method of calculating the Higgs mass is the best compromise between 
transparency and accuracy. 

Regarding the masses of the other four Higgses, the tree level results from [3] apply and are sufficient 
for the present purposes. We re-state them here. 

m 2 A o = 26/ sin(2/3) = 2y 2 + m 2 Hu + m 2 Hd (B19) 

m 2 H o = 7} ( m Ao + M| + ij (m 2 A o - M§ ) 2 + 4M|m^ 0 sin 2 (2fj)\ (B20) 

= + M\y. (B21) 

Appendix C: Application of the Checks and Iterative Procedure 

In this Appendix, we describe-for a single randomly generated initial point-two things: 1) the precise 
algorithm by which the checks described in Table III are applied and 2) the iterative numerical method used 
to solve for the B — L and SUSY scales. It is necessary to discuss these simultaneously since, as will 
become clear, they arc interrelated. We include this Appendix to give the reader insight into the details of 
our statistical method and to elucidate technical comments made in the main text. 

Before proceeding, it is helpful to note several things. A “point” here refers to a randomly generated 
choice of the parameters listed in Table II. For each point, we make working “guesses” of the initial values 
of AfsusY and Mb-l ■ These will be iteratively improved using a simple numerical method. For a fixed 
choice of randomly generated parameters and the two scales AfsusY and Mb~l specified, there is a unique 
solution for all of the RGEs and physical masses. That unique solution is found by our code using a 
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combination of analytic solutions (discussed throughout this paper) and numerical methods (not discussed 
in this paper). For the purposes of this Appendix, it is sufficient to know that the solution can indeed be 
calculated. It is also useful to note that, with the exception of the spill and convergence checks, the checks 
in Table III arc applied sequentially. For example, a point is subjected to the EW breaking check if and 
only if it passes the preceding B — L breaking and Zr bound checks. This means that a point that fails 
a particular' check a) has implicitly passed all previous checks and b) is immediately discarded and never 
subjected to subsequent checks. The sequential nature of these checks is what enables us to define the 
survival rates given in Table III. The spill checks and the convergence check, however, are different because 
they are not necessarily applied in a particular order and may even be applied multiple times to a single 
point. Nevertheless, if any point fails a spill or convergence check, at any step in the iterative process, 
we count that point as having passed all spill checks that appeal' above the failed check in Table III. This 
removes any ambiguity about how to define survival rates for the spill and convergence checks. 

Now we are prepared to discuss the main goals of this Appendix. For each randomly generated point, 
the initial guesses for Msusy and Mr-l are always taken to be 1 TeV and 2.5 TeV respectively. If the point 
with these initial guesses does not satisfy B — L breaking, then we count it as failing the B — L breaking 
check. If the point does not satisfy the Zr lower bound, then it is so counted. If it does not satisfy EW 
breaking, then it is so counted. If it does not satisfy the non-tachyonic stops check, it is so counted. 

If the guess for the B — L scale satisfies its definition, that is, if the RG calculation of Mz r (Mr-l) = 
Mr~l , to within 1%, and the guess for the SUSY scale satisfies its definition, that is, the RG calculation of 
m f 1 , rn t -2 satisfies y/m^ (Msusy ) m t 2 (Msusy) = Msusy to within 1%, then “convergence” has occurred 
and the steps in the next two paragraphs are skipped. 

If the guess for the B — L scale satisfies its definition to within 1%, then the rest of the steps in this 
paragraph are skipped. If not, the guess for the B — L scale is changed to Mz R . Using the same value for 
Msusy> and the new choice of Mr-l, we again run the RGEs for the same initial point. If Mz R not within 
1% of Mr-l, then the process is repeated. If the steps in this paragraph are repeated more than 300 times 
without success, then we count the point as having failed the convergence check. 

If the guess for the SUSY scale satisfies its definition to within 1%, the rest of the steps in this paragraph 
are skipped. If not, the guess for the SUSY scale is changed to y/m^ (Msusy ) m t 2 ( M$\jsy ) and we rerun 
the RGEs. If the point now does not satisfy EW breaking, it is counted as failing the EW breaking spill 
check. If the point now does not satisfy the non-tachyonic stops check, it is counted as failing the non- 
tachyonic stops spill check. If it does pass these checks, but Msusy does not satisfy its definition to within 
1%, then the steps in this paragraph are repeated. If they have been repeated more than 300 times without 


a conveniently chosen number which provides adequate opportunity for the iteration to converge. 



60 


success, the point is counted as failing the convergence check. 

Having successfully passed all of the previous criterion, we now must check the remaining checks. If 
the point does not satisfy the Zr bound, it is counted as failing the B — L bound spill check. If the point 
does not satisfy the sparticle bounds, it is so counted. If the point does not satisfy the Higgs mass check, it 
is so counted. If it does, however, satisfy all of these experimental checks, it is a valid point. 

The procedure described in the previous five paragraphs is represented pictorially by the “flow chart” in 
Fig. 18. 



FIG. 18. A “flow chart” showing how the checks are applied and how the iterative process of solving for the B — L 
and SUSY scales works. Every block that begins with the word “Check” has an outgoing red and green arrow. The 
green arrow is followed if the check is satisfied and the red arrow is followed if the check is not satisfied. 




































































61 


[1] S. Dimopoulos and H. Georgi, “Softly Broken Supersymmetry and SU(5),” Nucl. Phys. B 193, 150 (1981). 

[2] C. R. Nappi and B. A. Ovrut, “Supersymmetric Extension of the SU(3) x SU(2) x U(l) Model,” Phys. Lett. B 
113, 175 (1982). 

[3] S. P. Martin, “A Supersymmetry primer,” Adv. Ser. Direct. High Energy Phys. 21, 1 (2010) [hep-ph/9709356]. 

[4] C. S. Aulakh, A. Melfo, A. Rasin and G. Senjanovic, “Seesaw and supersymmetry or exact R-parity,” Phys. Lett. 
B 459, 557 (1999) [hep-ph/9902409], 

[5] C. S. Aulakh, B. Bajc, A. Melfo, A. Rasin and G. Senjanovic, “SO(IO) theory of R-parity and neutrino mass,” 
Nucl. Phys. B 597, 89 (2001) [hep-ph/0004031]. 

[6] K. S. Babu and R. N. Mohapatra, “Minimal Supersymmetric Left-Right Model,” Phys. Lett. B 668, 404 (2008) 
[arXiv:0807.0481 [hep-ph]]. 

[7] D. Feldman, P. Fileviez Perez and P. Nath, “R-parity Conservation via the Stueckelberg Mechanism: LHC and 
Dark Matter Signals,” JHEP 1201, 038 (2012) [arXiv: 1109.2901 [hep-ph]]. 

[8] P. Fileviez Perez and M. B. Wise, “Low Energy Supersymmetry with Baryon and Lepton Number Gauged,” 
Phys. Rev. D 84, 055015 (2011) [arXiv: 1105.3190 [hep-ph]]. 

[9] C. S. Aulakh and R. N. Mohapatra, “Neutrino as the Supersymmetric Partner of the Majoron,” Phys. Lett. B 
119, 136 (1982). 

[10] M. J. Hayashi and A. Murayama, “Radiative Breaking of SU{2)rXU{1)(b~l) Gauge Symmetry Induced by 
Broken N = 1 Supergravity in a Left-right Symmetric Model,” Phys. Lett. B 153, 251 (1985). 

[11] A. Masiero and J. W. F. Valle, “A Model for Spontaneous R Parity Breaking,” Phys. Lett. B 251, 273 (1990). 

[12] R. N. Mohapatra, “New Contributions to Neutrinoless Double beta Decay in Supersymmetric Theories,” Phys. 
Rev. D 34, 3457 (1986). 

[13] P. Fileviez Perez and S. Spinner, “Spontaneous R-Parity Breaking and Left-Right Symmetry,” Phys. Lett. B 673, 
251 (2009) [arXiv:0811.3424 [hep-ph]]. 

[14] V. Barger, P. Fileviez Perez and S. Spinner, “Minimal gauged U(1)(B-L) model with spontaneous R-parity 
violation,” Phys. Rev. Lett. 102, 181802 (2009) [arXiv:0812.3661 [hep-ph]]. 

[15] L. L. Everett, P. Fileviez Perez and S. Spinner, “The Right Side of Tev Scale Spontaneous R-Parity Violation,” 
Phys. Rev. D 80, 055007 (2009) [arXiv:0906.4095 [hep-ph]]. 

[16] R. N. Mohapatra, “Mechanism for Understanding Small Neutrino Mass in Superstring Theories,” Phys. Rev. 
Lett. 56, 561 (1986). 

[17] P. Fileviez Perez and S. Spinner, “The Minimal Theory for R-parity Violation at the LHC,” JHEP 1204, 118 
(2012) [arXiv: 1201.5923 [hep-ph]]. 

[18] P. Fileviez Perez and S. Spinner, “Supersymmetry at the LHC and The Theory of R-parity,” Phys. Lett. B 728, 
489 (2014) [arXiv: 1308.0524 [hep-ph]]. 

[19] P. Fileviez Perez, “New Paradigm for Baryon and Lepton Number Violation,” arXiv: 1501.01886 [hep-ph]. 

[20] V. Braun, Y. H. He, B. A. Ovrut and T. Pantev, “A Heterotic standard model,” Phys. Lett. B 618, 252 (2005) 



62 


[hep-th/0501070]. 

[21] V. Braun, Y. H. He, B. A. Ovrut and T. Pantev, “The Exact MSSM spectrum from string theory,” JHEP 0605, 
043 (2006) [hep-th/0512177], 

[22] V. Braun, Y. H. He and B. A. Ovrut, “Supersymmetric Hidden Sectors for Heterotic Standard Models,” JHEP 
1309, 008 (2013) [arXiv: 1301.6767 [hep-th]]. 

[23] A. Lukas, B. A. Ovrut, K. S. Stelle and D. Waldram, “The Universe as a domain wall,” Phys. Rev. D 59, 086001 
(1999) [hep-th/9803235]. 

[24] A. Lukas, B. A. Ovrut, K. S. Stelle and D. Waldram, “Heterotic M theory in five-dimensions,” Nucl. Phys. B 
552, 246 (1999) [hep-th/9806051]. 

[25] V. Braun, B. A. Ovrut, T. Pantev and R. Reinbacher, “Elliptic Calabi-Yau threefolds with Z(3) x Z(3) Wilson 
lines,” JHEP 0412, 062 (2004) [hep-th/0410055], 

[26] V. Braun, Y. H. He, B. A. Ovrut and T. Pantev, “Vector bundle extensions, sheaf cohomology, and the heterotic 
standard model,” Adv. Theor. Math. Phys. 10, 525 (2006) [hep-th/0505041], 

[27] M. Ambroso and B. Ovrut, “The B-L/Electroweak Hierarchy in Heterotic String and M-Theory,” JHEP 0910, 
011 (2009) [arXiv:0904.4509 [hep-th]]. 

[28] M. Ambroso and B. A. Ovrut, “The B-L/Electroweak Hierarchy in Smooth Heterotic Compactifications,” Int. J. 
Mod. Phys. A 25, 2631 (2010) [arXiv:0910.1129 [hep-th]]. 

[29] M. Ambroso and B. A. Ovrut, Int. J. Mod. Phys. A 26, 1569 (2011) [arXiv: 1005.5392 [hep-th]]. 

[30] B. A. Ovrut, A. Purves and S. Spinner, “Wilson Lines and a Canonical Basis of SU(4) Heterotic Standard 
Models,” JHEP 1211, 026 (2012) [arXiv: 1203.1325 [hep-th]]. 

[31] Z. Marshall, B. A. Ovrut, A. Purves and S. Spinner, “Spontaneous /AParity Breaking, Stop LSP Decays and the 
Neutrino Mass Hierarchy,” Phys. Lett. B 732, 325 (2014) [arXiv:1401.7989 [hep-ph]]. 

[32] Z. Marshall, B. A. Ovrut, A. Purves and S. Spinner, “LSP Squark Decays at the LHC and the Neutrino Mass 
Hierarchy,” Phys. Rev. D 90, 015034 (2014) [arXiv: 1402.5434 [hep-ph]]. 

[33] B. A. Ovrut, A. Purves and S. Spinner, “A Statistical Analysis of the Minimal SUSY B-L Theory,” 
arXiv: 1412.6103 [hep-ph], 

[34] D. K. Ghosh, G. Senjanovic and Y. Zhang, “Naturally Light Sterile Neutrinos from Theory of R-parity,” Phys. 
Lett. B 698, 420 (2011) [arXiv: 1010.3968 [hep-ph]]. 

[35] V. Barger, P. Lileviez Perez and S. Spinner, “Three Layers of Neutrinos,” Phys. Lett. B 696, 509 (2011) 
[arXiv: 1010.4023 [hep-ph]]. 

[36] B. Mukhopadhyaya, S. Roy and L. Vissani, “Correlation between neutrino oscillations and collider signals of 
supersymmetry in an R-parity violating model,” Phys. Lett. B 443, 191 (1998) [hep-ph/9808265], 

[37] E. J. Chun and J. S. Lee, “Implication of Super-Kamiokande data on R-parity violation,” Phys. Rev. D 60, 
075006 (1999) [hep-ph/9811201], 

[38] E. J. Chun and S. K. Kang, “One loop corrected neutrino masses and mixing in supersymmetric standard model 
without R-parity,” Phys. Rev. D 61, 075012 (2000) [hep-ph/9909429]. 

[39] M. Hirsch, M. A. Diaz, W. Porod, J. C. Romao and J. W. L. Valle, “Neutrino masses and mixings from super- 



63 


symmetry with bilinear R parity violation: A Theory for solar and atmospheric neutrino oscillations,” Phys. Rev. 
D 62, 113008 (2000) [Erratum-ibid. D 65, 119901 (2002)] [hep-ph/0004115]. 

[40] G. Gamberini, G. Ridolfi and F. Zwirner, “On Radiative Gauge Symmetry Breaking in the Minimal Supersym¬ 
metric Model,” Nucl. Phys. B 331, 331 (1990). 

[41] K.A. Olive et al. (Particle Data Group), Chin. Phys. C, 38, 090001 (2014). 

[42] S. P. Martin and M. T. Vaughn, “Two loop renormalization group equations for soft supersymmetry breaking 
couplings,” Phys. Rev. D 50, 2282 (1994) [Erratum-ibid. D 78, 039903 (2008)] [hep-ph/9311340]. 

[43] J. R. Ellis, K. Enqvist, D. V. Nanopoulos and K. Tamvakis, “Gaugino Masses and Grand Unification,” Phys. 
Lett. B 155, 381 (1985). 

[44] K. Choi, A. Falkowski, H. P. Nilles and M. Olechowski, “Soft supersymmetry breaking in KKLT flux compact- 
ification,” Nucl. Phys. B 718, 113 (2005) [hep-th/0503216]. 

[45] A. Djouadi, “The Anatomy of electro-weak symmetry breaking. I: The Higgs boson in the standard model,” 
Phys. Rept. 457, 1 (2008) [hep-ph/0503172], 

[46] N. Arkani-Hamed and S. Dimopoulos, “Supersymmetric unification without low energy supersymmetry and 
signatures for fine-tuning at the LHC,” JHEP 0506, 073 (2005) [hep-th/0405159]. 

[47] M. E. Cabrera, J. A. Casas and A. Delgado, “Upper Bounds on Superpartner Masses from Upper Bounds on the 
Higgs Boson Mass,” Phys. Rev. Lett. 108, 021802 (2012) [arXiv:l 108.3867 [hep-ph]]. 

[48] G. F. Giudice and A. Strumia, “Probing High-Scale and Split Supersymmetry with Higgs Mass Measurements,” 
Nucl. Phys. B 858, 63 (2012) [arXiv:l 108.6077 [hep-ph]]. 

[49] LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/04-01.1, 
http://lepsusy.web.cern.ch/lepsusy. 

[50] LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/02-09.2, 
http://lepsusy.web.cern.ch/lepsusy. 

[51] LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/02-10.1, 
http://lepsusy.web.cern.ch/lepsusy. 

[52] LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/01-03.1, 
http://lepsusy.web.cern.ch/lepsusy. 

[53] LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/02-04.1, 
http://lepsusy.web.cern.ch/lepsusy. 

[54] CMS Collaboration [CMS Collaboration], “Search for supersymmetry in hadronic final states using MT2 with 
the CMS detector at sqrt(s) = 8 TeV,” CMS-PAS-SUS-13-019. 

[55] G. Aad el al. [ATLAS Collaboration], “Search for squarks and gluinos with the ATLAS detector in final states 
with jets and missing transverse momentum using s/s = 8 TeV proton-proton collision data,” JHEP 1409, 176 
(2014) [arXiv: 1405.7875 [hep-ex]]. 

[56] [ATLAS Collaboration], “Search for high-mass dilepton resonances in 20 /6 _1 of pp collisions at sjs = 8 TeV 
with the ATLAS experiment,” ATLAS-CONF-2013-017. 

[57] CMS Collaboration [CMS Collaboration], “Search for Resonances in the Dilepton Mass Distribution in pp 



64 


Collisions at sqrt(s) = 8 TeV,” CMS-PAS-EXO-12-061. 

[58] G. Aad et al. [ATLAS Collaboration], “Measurement of the Higgs boson mass from the H —► 77 and H — » 
ZZ* —» 4£ channels with the ATLAS detector using 25 tb 1 of pp collision data,” Phys. Rev. D 90, 052004 
(2014) [arXiv: 1406.3827 [hep-ex]]. 

[59] S. Chatrchyan et al. [CMS Collaboration], “Observation of a new boson with mass near 125 GeV in pp collisions 
at yfs = 7 and 8 TeV,” JHEP 1306, 081 (2013) [arXiv: 1303.4571 [hep-ex]]. 

[60] J. R. Ellis and D. V. Nanopoulos, “Flavor Changing Neutral Interactions in Broken Supersymmetric Theories,” 
Phys. Lett. B 110, 44 (1982). 

[61] R. Barbieri and R. Gatto, “Conservation Laws for Neutral Currents in Spontaneously Broken Supersymmetric 
Theories,” Phys. Lett. B 110, 211 (1982). 

[62] B. A. Campbell, “Supersymmetry And Neutral Flavor Nonconservation,” Phys. Rev. D 28, 209 (1983). 

[63] J. R. Ellis, S. Ferrara and D. V. Nanopoulos, “CP Violation and Supersymmetry,” Phys. Lett. B 114, 231 (1982). 

[64] W. Buchmuller and D. Wyler, “CP Violation and R Invariance in Supersymmetric Models of Strong and Elec- 
troweak Interactions,” Phys. Lett. B 121, 321 (1983). 

[65] J. Polchinski and M. B. Wise, “The Electric Dipole Moment of the Neutron in Low-Energy Supergravity,” Phys. 
Lett. B 125, 393 (1983). 

[66] F. del Aguila, M. B. Gavela, J. A. Grifols and A. Mendez, “Specifically Supersymmetric Contribution to Electric 
Dipole Moments,” Phys. Lett. B 126, 71 (1983) [Erratum-ibid. B 129, 473 (1983)]. 

[67] D. V. Nanopoulos and M. Srednicki, “The Demon of Local SUSY,” Phys. Lett. B 128, 61 (1983). 

[68] L. J. Hall, R. Rattazzi and U. Sarid, “The Top quark mass in supersymmetric SO(10) unification,” Phys. Rev. D 
50, 7048 (1994) [hep-ph/9306309, hep-ph/9306309], 

[69] M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, “Electroweak symmetry breaking and bottom 
- top Yukawa unification,” Nucl. Phys. B 426, 269 (1994) [hep-ph/9402253], 

[70] N. Arkani-Hamed, H. C. Cheng and L. I. Hall, “A New supersymmetric framework for fermion masses,” Nucl. 
Phys. B 472, 95 (1996) [hep-ph/9512302], 

[71] C. F. Berger, J. S. Gainer, J. L. Hewett and T. G. Rizzo, “Supersymmetry Without Prejudice,” JHEP 0902, 023 
(2009) [arXiv:0812.0980 [hep-ph]]. 

[72] J. A. Conley, J. S. Gainer, J. L. Hewett, M. P. Le and T. G. Rizzo, “Supersymmetry Without Prejudice at the 
LHC,” Eur. Phys. J. C 71, 1697 (2011) [arXiv: 1009.2539 [hep-ph]]. 

[73] L. E. Ibanez and G. G. Ross, “SU(2)-L x U(l) Symmetry Breaking as a Radiative Effect of Supersymmetry 
Breaking in Guts,” Phys. Lett. B 110, 215 (1982). 

[74] J. R. Ellis, D. V. Nanopoulos and K. Tamvakis, “Grand Unification in Simple Supergravity,” Phys. Lett. B 121, 
123 (1983). 

[75] L. Alvarez-Gaume, J. Polchinski and M. B. Wise, “Minimal Low-Energy Supergravity,” Nucl. Phys. B 221, 495 
(1983). 

[76] W. Porod, M. Hirsch, J. Romao and J. W. F. Valle, “Testing neutrino mixing at future collider experiments,” 
Phys. Rev. D 63, 115004 (2001) [hep-ph/0011248], 



65 


[77] M. Hirsch and W. Porod, “Neutrino properties and the decay of the lightest supersymmetric particle,” Phys. Rev. 
D 68, 115007 (2003) [hep-ph/0307364], 

[78] P W. Graham, D. E. Kaplan, S. Rajendran and P. Saraswat, “Displaced Supersymmetry,” JHEP1207, 149 (2012) 
[arXiv: 1204.6038 [hep-ph]]. 

[79] P W. Graham, S. Rajendran and P. Saraswat, “Supersymmetric crevices: Missing signatures of R -parity viola¬ 
tion at the LHC,” Phys. Rev. D 90, no. 7, 075005 (2014) [arXiv: 1403.7197 [hep-ph]]. 

[80] CMS Collaboration [CMS Collaboration], “Search for Pair-production of First Generation Scalar Leptoquarks 
in pp Collisions at sqrt s = 8 TeV,” CMS-PAS-EXO-12-041. 

[81] E. J. Chun, S. Jung, H. M. Lee and S. C. Park, “Stop and Sbottom LSP with R-parity Violation,” Phys. Rev. D 
90, no. 11, 115023 (2014) [arXiv: 1408.4508 [hep-ph]]. 

[82] F. S. Queiroz, K. Sinha and A. Strumia, “Leptoquarks, Dark Matter, and Anomalous LHC Events,” Phys. Rev. 
D 91, no. 3, 035006 (2015) [arXiv: 1409.6301 [hep-ph]]. 

[83] B. Allanach, A. Alves, F. S. Queiroz, K. Sinha and A. Strumia, “Interpreting the CMS £ + £~jjpr Excess with a 
Leptoquark Model,” arXiv: 1501.03494 [hep-ph], 

[84] J. A. Evans and Y. Kats, “LHC Coverage of RPV MSSM with Light Stops,” JHEP 1304, 028 (2013) 
[arXiv: 1209.0764 [hep-ph]]. 

[85] D. Aristizabal Sierra, M. Hirsch and W. Porod, “R-parity violating sneutrino decays,” JHEP 0509, 033 (2005) 
[hep-ph/0409241 ]. 

[86] D. Aristizabal Sierra, D. Restrepo and S. Spinner, “LSP sneutrino decays into heavy standard model pairs,” 
JHEP 1305, 046 (2013) [arXiv: 1212.3310 [hep-ph]]. 

[87] A. Bartl, M. Hirsch, T. Kernreiter, W. Porod and J. W. F. Valle, “Testing the mechanism of R parity breaking 
with slepton LSP decays,” JHEP 0311, 005 (2003) [hep-ph/0306071], 

[88] M. A. Luty, “2004 TASI lectures on supersymmetry breaking,” hep-th/0509029. 

[89] I. Antoniadis, E. M. Babalic and D. M. Ghilencea, “Naturalness in low-scale SUSY models and ”non-linear” 
MSSM,” Eur. Phys. J. C 74, no. 9, 3050 (2014) [arXiv: 1405.4314 [hep-ph]]. 

[90] P. Ciafaloni and A. Strumia, “Naturalness upper bounds on gauge mediated soft terms,” Nucl. Phys. B 494, 41 
(1997) [hep-ph/9611204], 

[91] B. de Carlos and J. A. Casas, “One loop analysis of the electroweak breaking in supersymmetric models and the 
fine tuning problem,” Phys. Lett. B 309, 320 (1993) [hep-ph/9303291], 

[92] J. A. Casas, J. R. Espinosa and I. Hidalgo, “The MSSM fine tuning problem: A Way out,” JHEP 0401, 008 
(2004) [hep-ph/0310137], 

[93] J. R. Ellis, K. Enqvist, D. V. Nanopoulos and F. Zwirner, “Observables in Low-Energy Superstring Models,” 
Mod. Phys. Lett. A 1, 57 (1986). 

[94] R. Barbieri and G. F. Giudice, “Upper Bounds on Supersymmetric Particle Masses,” Nucl. Phys. B 306, 63 
(1988). 

[95] G. W. Anderson and D. J. Castano, “Measures of fine tuning,” Phys. Lett. B 347, 300 (1995) [hep-ph/9409419], 

[96] M. E. Machacek and M. T. Vaughn, “Two Loop Renormalization Group Equations in a General Quantum Field 



66 


Theory. 2. Yukawa Couplings,” Nucl. Phys. B 236, 221 (1984). 



